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FOREWORD 


This report, prepared by the Dynamics and Loads Section, Martin Marietta 
Corporation, Denver Division, under Contract NAS5-11996, presents the re- 
sults of a study whose purpose was to develop a computer program system 
for dynamic simulation and stability analysis of passive and actively 
controlled spacecraft. The study was performed from May 1973 to April 
1975 and was administered by the National Aeronautics and Space Admin- 
istration, Goddard Space Flight Center, Greenbelt, Maryland, under the 
direction of Mr. Joseph P. Young. 


The report is published in four volumes: 


Volume I 
Volume II 
Volume III 
Volume IV 


Theory 

Program Users 1 Guide 
Demonstration Problems 
Program Listing 
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ABSTRACT 


A theoretical development and associated digital computer program 
system for the dynamic simulation and stability analysis of pas- 
sive and actively controlled spacecraft is presented. The dynamic 
system (spacecraft) is modeled as an assembly of rigid and/or 
flexible bodies not necessarily in a topological tree configura- 
tion. The computer program system may be used to investigate 
total system dynamic characteristics including interaction effects 
between rigid and/or flexible bodies, control systems, and a 
wide range of environmental loadings. Additionally, the program 
system may be used for design of attitude control systems and for 
evaluation of total dynamic system performance including time do- 
main response and frequency domain stability analyses. 

Volume I presents the theoretical developments including a des- 
cription of the physical system, the equations of dynamic equi- 
librium, discussion of kinematics and system topology, a complete 
treatment of momentum wheel coupling, and a discussion of gravity 
gradient and environmental effects. 

The development of synthesis and analysis techniques for the 
linearized system includes a discussion of the numerical linear- 
ization technique, procedures for definition of system transfer 
functions, and linear time domain response. 

Volume II is a program users* guide and includes a description of 
the overall digital program code, individual subroutines and a 
description of required program input and generated program out- 
put. 

Volume III presents the results of selected demonstration prob- 
lems that illustrate all program system capabilities. 

Volume IV contains a listing of the digital code. 
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PROGRAM SYSTEM OVERVIEW 


This volume is intended to provide the reader with sufficient 
understanding of program system DISCOS*and its capabilities so 
as to permit a user to employ the program as a basic tool to 
analyze the behavior of a wide range of dynamical problems. 
Specific emphasis will be on a simulation for multiply - inter- 
connected spinning elastic bodies responding under the combined 
influences of external environments and either active or pas- 
sive control. 


A . INTRODUCTION 


The simulation employs a state-space approach that was developed 
in detail in Volume I.. The state-space formulation provides an 
attractive basis for simulation of nonlinear dynamical problems 
in a general sense as well as permitting linearization of the 
governing equations to provide an additional foundation with 
which to evaluate frequency domain and linearized time domain 
charac teristics . 

An attempt has been made to relieve the user from the require- 
ment of having to communicate with the digital program via large 
amounts of bulk data input. Although the program has many op- 
tions available, the program data stream has been organized to 
require only a minimal amount of basic input data for a particu- 
lar simulation. The data requirements have been further consol- 
idated in a manner that is quite definitive for the physical 
system being simulated. In summary, the user can quite easily 
relate to the particular elements of the program requirements 
and thus minimize setup time required to prepare data input 
for a given problem. In addition, a set of self-checking fea- 
tures has been included in an attempt to identify and check 
certain compa tabi li ties that are necessary for a proper simu- 
lation of a physically realizable system. 

In an overall sense, the digital program can be employed by the 
user to obtain 

1. nonlinear time response, 

2. interaction constraint forces, 


Dynamic Interaction Simulation of Controls and Structure 
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3. total system resonance properties, 

4. frequency domain response and stability information, 

5. linearized time response. 

The program outputs consist of printed and plotted results de- 
pic ting 

1. dynamic model construction, 

2. time domain response, 

3. frequency domain characteristics. 

The printed outputs are of a fixed form while the user controls 
the plotted information through the input data stream. 

B. SIMULATION OVERVIEW AND NOMENCLATURE 


This discussion identifies the basic nomenclature used to syn- 
thesize a typical assembly of interconnected bodies. The 
theoretical development, program users f manual, and demon- 
stration problems make extensive reference to various termin- 
ologies that are clarified here. Figure I.B-1 provides a visual 
display that illustrates many of the items being discussed and 
will be repeatedly referred to in the ensuing discussions. 

The overall system "topology" is identified by the user via the 
input integer array ITOPOL that contains the necessary informa- 
tion describing which "hinges" interface which bodies. Each 
body contains a body reference point that is the origin of an 
orthogonal cartesian body axis system. This point need not 
coincide with the body center of mass. 

Contiguous "bodies" are interfaced through a "hinge". We say 
interfaced in lieu of connected to emphasize the fact that the 
common "hinge" point between contiguous bodies may actually 
permit relative transla tional motion of the two bodies at the 
hinge. The degree of fixity at the hinge is identified by the 
user via the input integer array IHDATA. A typical body may 
contain "sensor" points that identify particular points where 
additional information is required to complete the desired 
simulation, A sensor point might sense on position or rate 
for the control system inputs, but could also represent a point 
on a body where certain other information is desired, such as a 
momentum wheel location or a point of force/torque application. 
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The integer input array IFTSMW identifies the body, where a 
particular sensor point is located. 


A body may also contain "momentum wheels". This special con- 
sideration accommodates a disk or rotating mass with a single 
relative rotational degree of freedom into the simulation with- 
out introducing another body. The momentum wheel capability is 
more efficient for the simulation of a single degree of freedom 
rotating mass than is constraining 5 of 6 rigid body degrees of 
freedom via constraint equations. All momentum wheels must 
have an associated sensor point. A wheel may either be active 
or constant speed; an active wheel has a variable spin rate and 
receives an input torque (generally via some sensor output re- 
lationship) and a shaft torque is applied to the wheel inducing 
a wheel angular acceleration. The array IMO identifies whether 
or not the wheel is active, and which axis is the spin axis. 

The reference axis for the wheel is the same as the sensor 
point axis system where the wheel is located. The array AMO 
identifies the wheel spin rates (initial rate only for the ac- 
tive wheel) and the wheel spin inertia about the spin axis. 

The system state vector is arranged in a specific manner within 
the program and it is necessary for the user to be very familiar 
with this arrangement for a number of reasons. First, the user 
must know where certain variables are located so that he can 
couple the control law into the simulation, and secondly, the 
user must know the order of the state variables in order to 
interpret results. Figure I.B-1 presents the state variable 
order consistent with the illustrative problem and other re- 
lated information. The state variables shown do indeed repre- 
sent a typical arrangement in that all of the various types of 
variables resulting 'from the multiple options available within 
the simulation are present. The order of the constraints (A) 
is also noted. Note that the user introduces the control varia- 
bles into the state vector but these variables (8) will always 
appear after the betas ( 13 ). Furthermore, the user may also in- 
troduce auxiliary variables (plant sensor signals and control 
system outputs) for use in the linearized studies. These auxil- 
iary variables should be placed (by the user) after the control 
variables (8) and in the order: plant sensor signals (Xg S ) 

followed by the control system outputs (B) . 



GENERAL USAGE INFORMATION 


There are a number of guidelines that must be adhered to in 

setting up a particular simulation. Some were detailed pre- 
viously and are concisely summarized here: 

1. there must be at least two bodies; a single body problem 
is simulated by including a dummy body that is not con- 
nected to the body to be analyzed; 

2. body no,.l is always positioned relative to the-inertial 
reference ; 

3. bodies are numbered from 1 to NB in an arbitrary order; 

4. every body (except body 1) must have at least one hinge; 
body 1 must have at least two hinges; 

5. hinges are numbered from 1 to NH in an arbitrary order but 
hinge no, 1 is, by definition, the hinge on body 1 between 
body 1 and the inertial origin; hence, hinge no, 1 can only 
appear on body 1; 

6. there must be a t least one sensor point for a given simula- 
tion ; 

7. sensor points are numbered from 1 to NS in an arbitrary 
manner; 

8. a typical flexible body requires mass and modal data that 
reflects a coordinate system that, is consistent with the 
body axis reference system for that body, e,g,, a modal 
coupling approach establishing modal properties for a given 
body would have to use the same reference body axis system; 

9. for frequency domain studies, there can only be as many con- 
trol output variables identified to introduce into the state 
equations as there are control system variables to begin 
with. Similarly, there can be no more sensor signal varia- 
bles identified than plant variables which appear in the 
original independent state equations. 

10. the user must make certain that the user supplied package 
has dimensions consistent with NHMAX for the arrays 



SK(NSK, NHMAX), DK(NDK, NHMAX), and HNGT (NHT , NHMAX) 
whore NHMAX = dimensioned maximum number of hinges, 

NSK = 3 or 6 depending upon nature of hinge free- 
dom, 

NDK = 3 or 6 depending upon nature of hinge free- 
dom, 

NHT = 3 or 6 depending upon nature of hinge free- 
dom, 

and if rotation only, then NSK = NDK = NHT = 3, if rotation 
and translation, then NSK = NDK — NHT = 6; 

11, the inertial properties of all the momentum wheels in a 
particular body must be included in that body s inertia 
description (whether rigid or flexible) since inertial 
coup ling is used . 




PROGRAM SEGMENTATION 


The digital code has been segmented into an executive overlay 
which governs the succeeding program flow and four supporting 
primary overlays, each with a separate and dedicated purpose. 
The basic program flow is depicted in Figure II.A-1. 



Figure II.A-1. DISCOS Program Segmentation 




Table II. A- 1 summarizes the intended purpose of the fundamental 
components in the program structure. 


Table II. A- 1 DISCOS Overlay Description 


1 

Primary Overlays 1 

Segment Name 

Purpose 

MAIN 

Executive program control 

DYNS10 

Data input 

1 

DYNS20 

Simulation of problem, lineari- 
zation of state eq’s, nonlinear 
time response 

DYNS30 

Plot results from nonlinear or 
linearized time response 

DYNS40 

Frequency domain analysis, linear- 
ized time response, frequency domain 
displays, (Bode, Nichols, Nyquist, 
Root Locus) 

Secondary Overlays (called from DYNS10) 

„ 

MSMODL 

Flexible body data inputs for j 

lumped mass representation j 

MSMODC 

j 

Flexible body data inputs for j 

consistent mass representation j 

j 


The executive overlay (MAIN) initiates the simulation by read- 
ing job iden tif ica ti on information and then passes control to 
the first primary overlay (DNYS10) which represents the basic 
data input segment. This overlay may be viewed as the program 
segment which builds the model from the input data. A series 
of topology checks are made as the data is loaded within this 
overlay to better assure proper modeling of the physical system. 
This overlay utilizes two additional secondary overlays for pro- 
cessing certain types of inertial and modal data. 



After overlay DYNS10 has structured the basic data for simula- 
tion, control is returned to the executive overlay which in 
turn passes control on to the second overlay DYNS20 . This 
overlay performs the actual problem mechanization and develops 
the nonlinear formulation which is the foundation for the en- 
tire dynamic simulation program. 

During a given simulation, the executive overlay always calls 
both the first and second primary overlays (DYNS10 and DYNS20) 
but, depending upon certain input control parameters, may or 
may not call the time history plot overlay DYNS30 or the lin- 
earized system analysis overlay DYNS40 . 

Simulation of a particular problem has its basis within the 
algorithms contained in the program subroutine YD0T which es- 
tablishes the canonical first-order differential equations that 
govern the dynamical motion. This routine in turn addresses 
another subprogram T0RQUE which in turn activates the user 
supplied modules that relate to the particular simulation being 
considered. These modules are discussed in further detail in a 
following section . 




III. 


DELINEATION OF USER- SUPPLIED MODULES 


The program has been written under the assumption that certain 
user- supp lied modules are available to complete a given pro em. 

In this manner, the user has considerable latitude with regar 
to how certain particulars related to a given simulation are to 
be handled. Control law specification, external torque inputs, and 
identification of plant sensor signals and control system out- 
puts are examples of items handled by the user. With this con- 
cept in mind, several subprograms have been placed under user 
control but with certain restrictions and guidelines to which 
the user must adher. Later comments will identify certain re- 
quirements associated with these user supplied modules. 


LOGIC FLOW 


It is worthwhile to consider a flow chart segment of the pro- 
gram (Figure III.A-1) and its chronology within the solution 
process. The order which the user supplied modules are called 
is indicated by the integers 1 through 7 (for subroutines) and 
8 and 9 for functions. 
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User-Supplied Subroutines 




Figure III.A-1 Chronology of Addressing User-Pak 
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B. 


SPECIFICS RELATED TO USER-SUPPLIED MODULES 


The separate user supplied subprograms each have. specific in- 
tended purposes and have been coded to fulfill these goals. 

The user can extend the scope of any of these modules with his 
own code, but there are certain items that these routines must 
perform. In any case, the potential user should be very famil- 
iar with many of the details of the user supplied package, and 
it 1 s with this fact* in mind that a separate discussion will 
now be devoted to c c i of the user supplied modules. Reference 
will be made to some of the programming logic contained in the 
DISCOS subroutine T0RQUE, and so this logic has been put into 
flow chart form as Figure III.B-1. Also, for reference pur- 
poses, the seven user-pak subroutines and two subfunctions 
(ADT, A DDT) are included in Appendix B corresponding to typical 
situations. 
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Creates external torques/forces 


Note: {G} is force/torque 

vector ^ RHS of equation 
of motion 

Adds momentum wheel shaft 
torques to {G} 

Adds hinge spring torques to |g| 
and sums spring energy to 
potential energy 

Gets AC due to thermal 
environment 


Now have RHS of equations 
of motion 


Creates additional equations 
for similarity transformation 
(used in linearization and 
stability package) 


Figure III.B-1 Subroutine T0RQUE Flow Diagram 
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1 . 


C0NTRL 


This is the first of the user-supplied routines and is always 
called by DYNS20 . The primary purpose of C0NTRL is tp esta- 
blish the time derivatives of the control system variables. 

These variables may be required by some of the other user rou- 
tines that arc activated after C0NTRL has been addressed. The 
routine must also establish the number of plant sensor signals 
(NXSS) and the number of control system outputs (NBTQ) which 
are transmitted through common block (/LDSIZE/) to the remainder 
of the program. For transfer function studies, the user is also 
required to identify whether or not transfer function poly- 
nominals are to be utilized. This is accomplished in a data 
statement (in C0NTRL) with the variable IsPLY which is the num- 
ber of polynominal ratio pairs (numerator and denominator) to 
be utilized. The first call to C0NTRL will read in the poly- 
nominal coefficients (for NPLY 4 0). 

Subroutine CONTRL contains a good deal of information pertain- 
ing to the simulation by virtue of its common blocks. Section 
C Identifies the constiuents of the common blocks contained in 
this and other modules. Additional common blocks can be esta- 
blished by the program user to transfer information between the 
separate user-pak modules. 

2. EXT0R 

This subroutine establishes the system external torques. 
Typically, this module can be utilized to accommodate such 
items as RCS (Reaction Control System) forces and torques, 
ae rodynamic s , and/or solar wind. The user can also extend 
this routine to include the addition of other state dependent 
torques. In summary, EXT0R, can be used as a M catch all for 
Inclusion of any additional forces and torques acting on the 
s vs tern. A single call to EXT0R from subroutine T0RQUE esta- 
blishes an integer array (ISNP) whose elements identify which 
sensor points are to be used for force/torque inputs. A vector 
containing torque and force components (ordered: T x , Ty, T z , 

F , F , , F 2 ) is then established for each of the force/torque 
sensor points and placed as a column into the array TEX. The 
vector of discrete forces and torques (referred to the local 
sensor-axis system) is returned to subroutine T0RQUE from 
EXT0R. These forces and torques are then transformed and added 
to the toLal system external force/ torque array ^g| . The user 
can bypass EXT0R related calculations by setting the variable 
NTEX equal to zero. 
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3. 


SHAFT! 


This routine establishes the shaft torque for each of the non- 
constant speed momentum wheels . Zeros are inserted for the 
torque contributions to the external torque vector for a con- 
stant speed wheel. 

4. KKLNGE 

This routine sets up hinge spring and dashpot torques/forces. 

It also accounts for potential energy contributions due to hinge 
spring deflections. The user must identify where spring rates 
and dashpot constants are to be found. This can easily be 
handled by a user specified equivalence statement within sub- 
routine KHENGE to locate the leading stiffness and damping 
elements within the data block identified as CNTDTA . Note : 
within subroutine KHINGE (see subroutine listing, Appendix B) 
there are statements of the following form 

DIMENSION SK(3 , NHMAX) , DK(3, NHMAX) , HNGT(3 , NHMAX) 


DO 10 1=1,3 
DO 15 1=1,3 


DO 20 


1=1,3 


where the integer 3 reflects the fact that consideration has 
been restricted to admitting only rotational springs at each 
hinge. If the user wants to also include springs/dashpots in 
relative translation at the hinge points, the three (3) in the 
statements above must be changed to a six (6), and appropriate 
spring rates and/or dashpots included within the data input 
array CNTDTA. Further, the equivalence statement locating the 
first spring rate ( SK( 1) ) and dashpot constant (DK(1)) reflect 
an order that is consistent with the hinge order; that is the 
first three elements in CNTDTA starting with the location cor- 
responding to the leading element of SK represents in order 
K#2 s for the first hinge. A similar relationship exists 
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5. 


for the array DK. 

Exam p le i In KHINGE note that 

EQUIVALENCE (CNTDTA(K) ,SK(1)) , (CNTDTA(L) ,DK(1)) 


and the array CNTDTA would by 


CNTDTA = 
element (K) 
element (L) 




^6 i ^9 2 ^^3 

Hinge 1 
(springs) 


^ei c e 2 c e : 

Hinge 1 
(dashpots) 


^ 9 ' * * 

Hinge NH 
(springs) 

c 8i c 9 2 C ^3 * ' " 

Hinge NH 
(dashpots) 


where the order of 9^ 2,3 *- s cons i- s t en t with the Euler rotation 
type (1-12) for the hinge q triad. 


The remainder of KHINGE is concerned with the proper placement 
of the spring/dashpot forces and torques onto the composite NB 
bodies (generalization of forces and torques) and should remain 
unchanged. 

The user can modify the referenced torques and forces immediate- 
ly after the (D0 10 L=1,NH) loop if he desires, but care must 
be taken to assure that the proper force or torque is correctly 
applied to accomplish the desired result. Several of the de- 
monstration problems (refer to Volume III) employ this process 
to apply control system outputs. 


GMISC 

This routine is reserved to implement torque/force contributions 
from thermal gradient effects. The entire state vector, along 
with component position and attitude information, is available 
via transfer through labeled common arrays. Section C provides 
more insight into the information contained in these common 
blocks. 


6. EQADD 

This routine establishes additional equations for use in the 
linearized time domain analyses. It must identify the number 
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of additional equations introduced via the variable NAUX (num- 
ber of auxiliary equations). These equations relate plant 
sensor signals, X^ s and control system output forces/ torques , 
B^, to the system state and in the specified order. The addi- 
tional variables must be placed in the state vector as the last 
NAUX state variables and in the order , then B*- and they 

become an integral part of system transfer function evaluations 

7 . LT0RQL 

k 

This routine establishes the portion of the right hand 

side of 


zi = A ij2 J 


+ b 



which is used for the linearized time response. This corres- 
ponds to the external excitations for the transformed variables 
zi, leading to evaluation of the perturbation response. 

8 . APT (Subfunction) 

This function is used in conjunction with ADDT to implement 
prescribed kinematical motion in the hinge coordinates. With 
reference to Figure I.B-1, the array IHDATA(I , J) ,I> 1, may have 
2 as an entry indicating that the Jth hinge has velocity and 
acceleration prescribed in that coordinate. The argument of 
this subfunction is: ADT(IC,T), with IC=6*(J-1)+(I-1) corres- 

ponding to IHDATA(I,J)=2. The integer IC and the time T are 
passed into ADT via argument by the calling subroutine so that 
the velocity (a) may be established for the proper hinge coor- 
dinate as a function of time. 

For a given rheonomic constraint, we note that both a and a 
must be set by subfunction. Now, it is conceivable that the 
user knows a as the exact mathematical time derivative of «. 

It would seem that the natural thing to do would be to create 
ADT and ADDT subfunctions to return consistent a and « respec- 
tively. This is not the best thing to do, however, because of 
numerical integration characteristics. The numerical integra- 
tion of ju} reflects the use of «. The resulting |u| reflects 
a numerically integrated « which cannot be consistent with a 
value obtained any way other than numerical integration. The 
consequences of this are seen as slight errors in motion res- 
ponse, but also, there is a large spurious change in system 
momenta . 
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The best wa y to effect rheonomic constraints is to use values 
of ft obtained from numerically integrating a. This can be 
easily done by using additional differential equations that 
are accommodated in the state vector as additional "control 
variables" or { 8 } . Thus, after all of the actual control varia- 
ble rates are established (in subroutine C0NTRL). one need only 
code additional expressions to set 8 (additiona 1) = “(desired). 
The statements within subfunction ADT merely return ADT=Y(K) ; 
the state vector Y is available in labeled common. /VECTOR/ , and 
K corresponds to the location in Y where the 5 = « control var- 
iable resides. Of course, IC must be tested such that the 
appropriate “ = 8 is returned. 

9 . ADDT (Subfunction) 

This function is discussed with regard to its relationship to 
ADT in Section (8) above. This function has arguments: 

ADDT (IC ,T) , exactly the same as ADT, and returns values of « 
for appropriate IC and T consistent with the “ returned by ADT. 
■Note in Figure III. A- 1, the chronology is such that subroutine 
C0NTRL is addressed prior to function ADDT. This is so that 
C0NTRL can establish a value of 8 (additional) = “(desired) 
to put fn the state vector time derivative (YDT, also available 
in labeled common/ VECTOR/) . Now, for the appropriate time T 
and IC, it is only necessary to set ADDT=YDT (K) , where again K 
corresponds to the location in Y where the 8 = “ auxiliary con- 
trol variable resides. 


C. DISCUSSION OF SELECTED COMMON BLOCK INFORMATION 


The program user will very often have a need to access certain 
information that is calculated and stored within the program in 
order to compute specific variables required for the user sup- 
plied modules. Such information about the simulation is stored 
in multi-dimensional array form within labeled common blocks. 
These data provide a good supplement to the state variable con- 
tent which has been previously discussed in that the user can 
extract both total and relative positions and rates for any 
component of the simulated dynamical system once he has a firm 
understanding of where certain data reside within the program. 
The following subsections will discuss selected common block 
arrays to better familiarize the potential user with their con- 
tent . 
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1. Common b lock/BHBSRD/ contains three separate groups of in- 
formation which the user may need to access. This information 
is concisely summarized in double and/or triple subscripted 
arrays as 

BS (6 , 6+NMDB0D,NSPMAX) 

R0L(3 , 3 ,NBMAX) 

D0L(3,NBMAX) 

where the following items are noted - 

NMDBOD = maximum dimensioned number of modes per body, 

NBMAX = maximum dimensioned number of bodies, 

NSPMAX = maximum dimensioned number of sensor points. 

The array, BS(i,j,k), contains the kinematical coefficients for 
all of the ''sensor" points. The rows (subscript i=l,2...6) of 
the array refer to (in order: u x , wy 3 a> z , u, v, w) the com- 

ponents of absolute angular and translationa 1 velocity (sensor 
referenced) at sensor point k. The columns of the array (sub- 
script j) refer to the j=l,2,...6 + no. of elastic modes on 
body containing sensor point k. Thus, in general, if we want 
to know the projection (the i th velocity component) onto the 
triad located at sensor point k, the following expression is 
noted 


Vel. = BS(i,j 1 



j 


The array, ROL(i,j,k) contains the rotation transformations 
relating the body axis systems to the inertial reference. The 
elements of the array are the direction cosines between the 
body axes, £j c , and the fixed inertial system, £ 0 . Subscript 
k denotes the body number. 

The array DOL(i,k) contains the three vector components, (X, 

Y, Z) , from the inertial reference to the body axis system, 
for each body. 

2. Common block /SPECIF/ contains information which the user 
may require. These arrays are 


BETAH(6 , NHMAX) 

BETAHD(6 , NHMAX) 

RS(3 , 3 , 2* (NSPMAX) ) 

DS (3 , 2* (NSPMAX) ) 

where the following items are noted - 

NHMAX = maximum dimensioned number of hinges, 

NSPMAX = maximum dimensioned number of sensor points. 

The arrays BETAH(i ,j ) and BETAHD(i.j) contain the hinge BETA's 
and rates respectively (for hinge j). The order (i subscript) 
is given as 



where is the i_th Euler angle rate consistent with I TYPE for 
hinge j and is the ith velocity component of point q relative 
to point p in the p frame for hinge j . 

The array RS(i,j,k) contains the rotation transformations 
(direction cosines) between the sensor point axis system and 
the body axis system (body on which sensor is located). wo 
sets of transformations ore identified for a given sensor 
point. The first represents misalignment of the two triads 
without elastic deformation and the second includes the elastic 
deformation. The ordering (subscript k) proceeds as follows: 
the /th sensor rotation (without elastic deformation) is located 
at k = 2*H~ 1 . The total rotation transformation for the /_th 
sensor is located at k= 2 *£ . For a rigid body, these two 
transformations are identical. 



The array DS(i,k) contains the three components of the vector 
from the body axis system to the body sensor points (in the 
body axis system). Here again, there are two sets of vectors 
(rigid body and rigid body + elastic) for each sensor point. 

The first is for rigid body and the second includes the elastic 
deformation. The addressing algorithm is the same as for 
RS(i,j,k). 
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IV. 


PROGRAM INPUTS 


The dynamic simulation program utilizes some basic data input 
5 17 ou t i ne s in an attempt to standardize a large amount of the 
bulk data input. Additional formatted inputs have been used 
where it is more meaningful (and more efficient) to do so. As 
will be noted in the following section, there is a large amount 
of data input via subroutines READ and READIM. Therefore, it 
is useful to familiarize the reader with these two routines 
prior to describing overall program data input requirements. 


A. DISCUSSION OF SUBROUTINES READ AND READIM 


These two subprograms are structured to read matrix arrays in 
floating point (real) notation (subroutine READ) or fixed point 
(integer) notation (subroutine READIM). A thorough discussion 
of the routines and their supporting subroutines is contained 
in Appendix A. The following discussion gives a cursory over- 
view of their usage . 

The routines arc activated by a FORTRAN call of the form: 


CALL READ (A, NR, NC , KR, KC) or 

CALL READIM (IA, NR, NC, KR, KC) 

where the arguments in the call statement are 

A, (IA) = floating (fixed) matrix array of size NR by NC 


NR 

NC 

KR 

KC 

A call to 
be in the 


number of rows in array 

number of columns in array 

row dimension of array in calling program 

column dimension of array in calling program 

either of these input routines requires that the data 
following format: 


IV -1 



1 . 


Subroutine READ 


First card - matrix name, NR, NC with format (A6,I4,I5) 
Middle cards - data with format (215, 4D17.8) 
first 15 is row number 

second 15 is column number of leading D17.8 field 
next 4D17.8 are elements of the array 

Last card - ten zeros in columns 1 through 10 

2 . Subroutine READIM 

First card - matrix name, NR, NC with format (A6,I4,I5) 
Middle cards - data with format (215, 1415) 
first 15 is row number 

second 15 is column number of leading 15 field 
next 1415 are elements of the array 

Last card - ten zeros in columns 1 through 10 
B. INPUT DATA STREAM 


This section presents the program system input data stream 
together with the data input control logic. The approach taken 
herein is to first introduce an overview of the data inputs and 
program control logic in the form of a flow diagram (Figure 
IV. B- 1) and- to then identify the details in much the same way 
as the FORTRAN code accepts the data inputs. This method of 
presentation has been chosen as it most closely relates to the 
actual processing of the user inputs for a given simulation. 

In addition, the user can follow the program control or switch- 
ing logic to determine just what data are required to complete 
a particular simulation. 
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Figure IV. B-l 

Program System DISCOS Data Stream FZow (Sheet 1 of 9) 
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Program DISCOS Data Stream Flow (Sheet 3 oj D) 
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Fijure IV. a- 1 

Program Sjstem DISCOS Data Stream Flow (Sheet 4 of 9) 
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Figure IV.B-1 

Prog mm System DISCOS Data Stream Flow (Sheet S of 9) 
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VECTOR ST7E 1 8Y 3 CONTAINING I NTF r - ER 
CONTROL DATA 

TODfTflfH = PR T NT CCNT°OL FOR time RE c p CNF E - P° T N'T 
EVERY IPDAT A 1 1 > MULTIPLE^ CE n Fl T P T 
TFCf*TA(2) = PLOT CONTROL FOR TIME PESFCNSf - SAVE 
EVERY TPOATA(2t MULTIPLE* CE CH.TAT 
TPDATBfT) = G PEP p C c ’ w NONLINEAR TIME RESPONSE 
- j p P R E CP M LINEAR TIME C® FREOUENOV 
FF S r C NiSE 


CALL READ (CNTOTAt It NONPAR, 1, YC.PNT) 

VECTCR SI2E 1 BY NCNPAR FCR CONTROL SYSTEM v/ARIA°LES 
ANO USFP SUPPLIED V A RI A°LFS 

THI R T S A CATCH-ALL VECTOR FO* 5 ti'E® 0 A K DATA 

INFORMATION IS PUT INTO CCMMCN /CCNTRL / ANT IS 
ir>ENTIFIE f ' R Y USE° SUPPLIED EQUIVALENCE Mfl R 
IN USER c A'*' 

FIRST NDFLTA ELEMENTS MUST CONTAIN INITIAL VALUES 
FOR C C LTA VARIABLES IN STATE VECTOR 

AOOITTCNAL SPACE TS AVAILABLE TO THE LS£R 


CALL READ CWV, It 4, 1, 5) 

VECTOR SI7E 1 BY 4 FOP GRAVITY GRADIENT OATfl 

WV(i> = PROJECTION CE GRAVITY VECTOR CN * INERTIAL AXIS 

WV ( 2 ) - PROJECTION OF GRAVITY VECTOR CN V TNE RT T At AXIS 

WV(3» = RPCJECTIOM DF GRAVITY VECTOR CN 7 INERTIAL AXTS 

WV(4> = RAOILS VECTOR ERCM GRAVITY SCU C CE TO GENERAL 

VICINITY OF BODY CLUSTER 


NOTE I F( w\M 11 «■ WV(21**2 ♦ WV(3)**2) 

EC Ot WVC4) IS IGNORED 
NF 3, WV<4> MUST PE GT 1 

NOTE — GRAVITY VECTOR COMPONENTS IN 
UNITS OF ACCELERATION 
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K'^UT FOP RIGID Q CDY 
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r* subroutine m^igid - * 
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p.***#**##***####****#******* ******************************************** 



CALL R c A r <V, 1, 4, 1, F) 

VECTOR SI 7 E 1 PY 4 

V ( 1 ) = MA C S CF PCOY 

VC?) = y POCY 9fF POINT TC °OOY CO, «CDY T°IAO) 

V f 3 > = Y (BCOY PEF POINT TC B^PY CO, P^OY Toifin) 

VfO = 7 ( n 0rY P£F POINT TO PP P Y PG, POOY "^lAO) 


■PALL C E AO (V, 1» 6, 1, 6) 
VECTCF BITE 1 °Y f 


* / xy dm, ect, 

•/ Tf- 1 


Vtl) - JYX ^Vol 

V( 2 ) = JYY 

V ( 3 ) = J 7 7 { c O OY PEFEPENCEO INEPTIAP - EO^Y tdtaO) 

V(4> = JYY (TNFPTTAS ABOUT PEF PCIM, NOT ®OOY CG) 

VC C ) = JX 7 
V( c > = JY 7 

jsjmo = NUYRFF CF MTNGFC CN POOY N - PYCLUSIVE CF HINGE 1, 0 C!)Y 1 


pr i: T=i,NH° 


0 . — ... FFAOtMT , FORMAT - 215) NO*-, ITYPE 
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c NCH = HINGE NUNOcP 

r ITYPE = EULER POTATION TYPE TO ORIENT HINGE 

c T Cl AC WPT F COY TRIAD 
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SUBROUTINE. MSMOOL - INPUT FOP FLEXTPLS FCQv, LUMPED MASS MATRIX 


-rALL READ (A, NJ, 1, KJCTNT, K^CC^) 

MATRIX SIZE NJ «Y 1 WHPPc nj 
CN ECCY N 

FOR TMC TTH JOINT - 
Afl,i» = JOINT MARR 


= NUMP£p OP JCTNTO 


■TAIL READ (A , N J , 6 , K JOT NT, K'lOCE) 

MATRIX SIZF NJ BY 6 

FC P TMP ITH JCTNT - 

A T 1 , 1 » = JOINT INERTTA, JVX 
A(T,2) = JOINT INERTIA, JYY 
ACT, 3) - JOINT INERTIA, J7 1 
A ( T » A ) = JOINT INERTIA, JXY 
A 1 1 , 5 1 = JOINT TN£RTifl t j X t 
A (1 , 6) = jn I NT INERTIA, JYZ 


J xy " / dm * ect * 
*vol 


On s k=i,2 


•CALI READ CA, NJ, 3, KJOTNT, KMODE) 

MATRIX SITE NJ qy 3 
FOR THE ITH JOINT - 

* = 1 A ( T , i I = JOINT STATIC MASS MpM£NT, SX 

A(l,2» = JOINT STATIC MASS MOMENT, c Y 

A ( I » 3 1 = JOINT STATIC MAPS MOMENT, SZ 


y-2 A ( 1 , 1 1 
A(I,2) 
MI , ?> 


= X ( 6 0 DV REF POINT TO JOINT, BODY TOlAf!) 
- Y (BOOY REF °OINT TO JOINT, Bony TPTAO) 
* Z (BOOY RPF POINT TO JOINT, pony TPT An ) 


5 CONTINUE 
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C CALL Pfflr (A, NJ , NE, KJOINT, KMCOF) 

r, 

C * A T B T. X ST7E KJ BY NE WHEPf NE = MUMPER OF c LASTTC, 

c irr£$ RETAINED ECO «OnY N 

C 

c FOE THE ITH JOINT - 

r 

r K-l A<I,J* - X niSPL ACt^ENT AT JCIM, M OCE J 

C k = 2 A(I,J> = Y n IS°l A CEMENT AT JOINT, MOTE J 

r K = 3 A 1 1 , J) = Z DISPLACEMENT AT JOINT, *COF J 

C i<=4 = THETA Y POTATION AT JCTNT, , 'CO r J 

r \f-5 A ( T , J) = THETA Y POTATION AT JOINT, m(30E J 

r K=E A ( I , J) = THETA 7 POTATION AT JOINT, MCCE J 
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C 10 CONTINUE 

C 
0 

c DO 2G K-l, 2 
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CALL PEAO (A, NE, NE, KJOTNT , YMCCE) 
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c MATPIX SI2E NE 6Y NE 
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r K"= 1 A - MODAL STIFFNESS 
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C K = 2 A = •'COAL CAHPINF 
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C 20 CONTINUE 
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c 

-OEACINIT, FORMAT * 8D1C.3) ( A I J» , J-i, NE) 

C 

c VEO TCP OF INITIAL MCDAL DEFLECTION COORDINATES 
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r RE AD < NIT, FORMAT = 8010. 3> t At J> , J= 1 ,N'E ) 

C 

c VFCTCP OF INITIAL MOOAL VELOCITY CCOPOTNATFS 

C 

c 
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FOLLOWING EULER ANGIES MEASURED 
IN UNDEFORHEO CONFIGURATION 


MHO = NUMBER CF HINGES CN POOY N - EXCLUSIVE CF HINGE 1, »OOY 1 
OC 150 L * 1 » N H P 


FE APCMT, F CR ► AT = 315) KPH, ITYPE, JC T N T 
NOH = HINGE NUMBER 

ITYFF * EULER PCTATICN TYPE TO ORIENT HI NGF 
TPIAC WRT PQOY TPIAO 

JCTNT = JCTNT NUMBER OCPRESRCNCING TC HINGE POINT 


P£AD(NIT, FORMAT - 3G10.3) ( WV I J) , J=l, 3) 
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= theta 1 

(FIRST 

ROTATION) 
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NS« = NUMBER 
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CN . BCOY N 
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IF «NSO ,EO. 

0) RETURN 




PH lft D L*l* NS° 


PEAOINIT, FORMAT = 115) NOS, ITYPE, JCTNT 
NCS s SENSOR POINT NUMPER 

ITYPE s EULER ROTATION TYPE TO ORIENT SFNSCR POINT 
TPIAO HPT FCOY TRIAD 

JOINT » JOINT NUMBER CORRESPONDING TO SENSOR POTNT 
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PEAOtNIT,FCPHAT = 3010. 3) (WV l J) , J=l, 3> 

EUIFP ANGLES TO CLIENT SENSCP POINT TPTAO 
0«CER OEFINEO PY TTVPE 

MV(1) * THETA 1 IFIPST DOTATION) 

WV ( 2) * THETA 2 (SECNfl ROTATION) 

WV (31 * THETA 3 THIRD ROTATION) 

160 CONTINUE 

RETURN 

ENO 


PEPHUTATTO 
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iniao = 3 TP MODAL DAMPING MAT°TX 



■n n vi i m i n o .t »-> n -j n o n o 


r 

r 


CALL READIM (JV, 1, NNCDT, 1, KAO) 

VF.CTOR SI7F 1 RY M * C H T = NUMPE C C p RICIO apCY NOCES * 

MJM D f° OF ELA^TTC MOFES 

JV(JJ = INC WH£?E TNO IS CCLU^K NO. WHTCH CFLfJ) OF C°TC 
MORAL “ATPI* will APPEAR t N 
r FVISEF M30E MATRIX 
TNC ST J REPLACE COLUMN 

INF c Q 0 P c ltTE COLUMN 

IMC LT 0 REPLACE COLUMN, CHANCE STPNC 


r 

r, 

r 

r 

r* 

C r 


r 

r* 

$ 

r 

o 

r 

r* 


c 

r 

r 

C 

n 

c 

r 

c 

r 

0 

r 

r 


pall 15 c A C (A, N 0 A , NO A , KA D , KAO ) 

**AT«;ix c T 7£ E* NX oy c,*NX CONTAINS CCN c T ST c mt 
|*C $S P FPC ESENT ATIOK 

°CV'-C,CLUKK COFFOlNATE CR0E9 MUST B r CCNSTSTr^r 
WITH THE C-:0S££ OF F9FE00M T A°l E , J n CF 

FALL p^AD (A, KOA, h M C 0 T , KA°, KA«) 

^ AT f I y FI7E g*NX ny n’mcot CONTAINS MFOAL 0 rr TNTT I ON 

THf CCIUMNP r^CFE CR^FP) WILL °F PFOorcrpcp qy 
The T h F IJ T >»PCTCF J \i 

TF ( I 0 1 AK , FC , : • AN D. T C I A F .EF. C) GC TC 11 

rtl L =FAO (F.H2, i, h^CDT, 1, KA PI 

VFCTCC SI2r i pv NHCCT CONTAINING SCLAVS (>f NATURAL 
FR r OU c NCI ^S 

0M2(J) = SC U A® E ff JTh NATU c AI FPEOUENFy f fOOESo r N'n I\'G 
TO JTH Incut mode SHAPE 

11 rCNTTNUE 


ORIGINAL PAGE IS 

op poor quality 
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r 

c 

c 

o 

c 

C I f (IFRR* . EG • C> GC TO 5 

C 

c 

r----- REAOTNIT , FORM AT = 16) JTYFCl 

p JTYPCL = RE F £Pf NO r JOINT NUMRER WHOSE GEO-FT°TC 

c POSITION 000° Cl NATES WILL °E U c E r To 

r ESTA°l ISH FIG 10 «00Y MCCAl * A T F IV 

C 

C «EA0(NIT,FCRKAT = 301C.7) ( OM? ( J ) f J=1 , 3 ) 

r ovgfn = x component nr vector tmt locates jtypol pt 

C OV2<2> = Y CC-PCNFNT OF VECTQP THAT LOCATE* JTYPCL p T 

n 0-2(3) = 7 COMPONENT OR VFOTO° THAT LOCATES JTYPCL PT 

C 

c VECTOR 00 -POKE NTS FEF£ c £0 TO POOY TPIAC 

C 

C 5 CONTINUE 

r 

r, 

c IF (ICIAK .EC. 1) GC TO EC 

C 

C CALL R£AD (A, NPA, NCA, KAO, KAP) 

c 

c STIFFNESS -ATPIX -- SEE NCTE RELOW 

r 

C 5 C CONTINUE 

n 

C IF (IOIAO .EG. 1> GC TC 6C 

C 

c 

C CALL READ (A, NP A, NCA, K A P , KAR) 

r 

o DAMPING NATPIX -- SEE NOTE REICH 

C 

C GO TO 61 

C 

0 63 CONTINUE 

C 

C 

C„ FEA0(NIT ,FCP YAT - flC10.3) (0M2 ( J) , J=1 ,NE) 

Q 

m VECTOR SI 7E 1 °Y NE = NUM°£R OF ELASTIC MOOES 

C RETAINED VIA INPUT JV SELECT TON VECTOR 

Q 

r. CN2(J) = NOOAL DAMPING RATIO F C° JTH ELASTIC MCOF 
C 

c .61 CONTINUE 

C 
C 
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KjrTE -- FOLLOWING EULER ANGLES M EASU^EO 

IK ukoefophed configuration 


C 

r 

r 

j 

c- 

c 

r 

C 

c 

c 

c 

c 

c* 

c 

0 
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r 

c 

c 

c 
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r 
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o 

c 

c 
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C 

c 

c 
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c 

0 
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r 

r 

c 

c 

c 


MH q = CP HINGES ON °OPY N - EXCLUSIVE CF W TNGF 1, POOV 1 

no lie 1 = 1 , NH0 


Ft ADCMTjFCFNAT = II?) NOH, TTYPE, JOINT 
NCH = HINOF. NUM“ER 

IIYPF. = EULER ROTATION TYFE TO ORIENT HTNGE 
TRIAO WRT ®ODY TRIAD 

JOINT - JOINT NUHRER OORRESRCNDTNn T r VT^ge °CTNT 


-..-0£fln(Kt T « FC°H AT = 3C10.7) ( OM? ( J) * J=1 » 3 ) 

EUIEF ANGLES TO OPIENT HINGE TRIAD - PERMUTATION 
CFCEF OFF INEC ev I TYPE 

= TMRTA 1 (FIPST ROTATION) 

0 " 2 < 2 > = THETA 2 (SECNO ROT ATT rw ) 

0*2 m = TH£ta 3 (THIRD RCT ATICN) 

110 rCNTIMJ r 


KP8 = NUHEEP CP SENSOR POINTS ON RCOY N 
TP ( NS° .EO. C) RETURN 
DO 12C L=1,NS P 


c £AO(NTT,FCRNAT = 315) NOS, I TYPE, JCTNT 
NCS - SENSOR POINT NUMREP 

ITYPE - EULER ROTATION TYPE TO ORIENT SENSOR POINT 
TRI AO WRT °0DY TRIAD 

JOINT = JCTNT NUHRER CORRESPONDING TO SENSOR POINT 
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MI SIN 3101 3d 300 d0i.tfa3wfN - ( 2 * I ) A 1d3 
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Ot Oa 00 (0 * 03 * AloN) al 

OVda 01 SAN313133 300 NOIiONOd 
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VIVO 0 31 IddOS assn V Ac OdMldBO SI A ldM 
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Oi UsO ASal3 3Hi NO OVda Si (ANV 31) VxVU 
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(ANV 3 1 ) VIVO X V d daSH IfldNl - 32SNAO 3Nlinad6nS *0 
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3 

0 


0 
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L. UOl-'O C- 


r w 

C« SUPRCUTINE OYNS30 - INPUT FOP TIME HISTORY PLOT OUTPUT 

C* 

C*********”**”«”*”*””*”**””"***** w ”***”**** ,w ”**” w * 

c 

c 

c PEAOfNlT, FORMAT = 1CA8) ( ICTITH I) , 1= 1 , 10) 

r. 

c BJ CHflcarTE» mastfr plot title will affeap 

c OK ALL FRAMES 

c 

c 

C PEAn(KIT,FCRHAT = 15) NSET 

r 

c NSET = NT. OF PICT SETS TO PE 

c OYCLEC THPOUGK EACH ^FT 

C IS ITMTTFr TC 16 'EFAFATE 

c VARIABLES TO **F. SFLEOTEO 

C FROM THE SET OF VARIABLES 

r; WCITTEN PE SU ° R CUT INF PLO' , WO 

C OF LPLTWP. 

0 

C 

C nn uso TSET = 1,KSET 

r 

C 

c — FFAOCMT ♦ F C R M. A T - 15) JPL 

C 

C JPL = NO, OF VAPIA°LES TO «E 

f* S r l EOT EC FPCM THE NCFI.CT 

r VARIABLES PRFVTCUSIY WPITTFN 

C ON NT A PE 3 FOR p L OTTTNG . 

r (JPl .LE. IF) FOP A SET 

C 

C NCFLOT = 1 + ?*NEC*NLAM *N U* JR ♦ 3 ♦ 3 *2 * N' Q + 5 

r (NONLINEAR ANALYSIS 

r, 

0 

r. N C FI OT = 1 + 2*NFC 

C (LINEAR ANALYSIS) 

r 

C WHERE 

0 

r. NFO = NO. C p EOLATIONS 

n TNTEGRATEO. 

C 

q NLA* = NO. OF LA m POA VARIABLES 

r 

C NU = NO. CP U VAPTA°LES 

c 

c JF = 6*NF 

r 

C NP = NO. CP OCCIFS. 

C 

r, 
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O O ‘-1 


NOTE 


MP 

* ** * * 


c 

r* 

c 

r. 

c 

c 

c 

c 


c 

c 

n 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

n 

c 

c 

c 

r 

c 

n 

c 

c 

c 

n 

c 

c 

c 

n 

c 

c 

n 


ne c 


_ * IPGFIX(J) ♦ NU ♦ NFFTA ♦ NHFlTfl 

* 

* ♦* V * 

J =1 


NU 


f'F 

- * IRGFLXl J) + 

* 

***** 

J=1 


KCF^r 

** *** 

* 

♦ iMCm 4- 6*N° 

* 


K =1 


KOFTfi s Sin* CF NUPRFR TF ZEROS ♦ ?UM rr 
NUMRFR of TWCS in ROWS 2 THPU 7 
CF ARRAY IHOATft 


Nl A" 


<5U* OF NUN°ER CF CNES ♦ SU* OF 
KUKPEF OF TWCS IN FCWS 2 THRU 7 
Or ARRAr IHQATA 
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.1 o o .t rt fi o o n o o o n o rt o >1 f) .1 n >i .1 o 


ORDER OF V AR I A°LES AND 
SIZE FC* A SINGLE PECOPO 
(FCP A SINGLE TINE, T> . 


NONLINEAR 

analysis' 


r 

r 

r 

r, 

c 

c 

c 

0 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r 

C 

c 

r 

r 

c 

c 

C LINEAR 

C ANALYSIS 

r, 

o 

c 

o 

c 


V AR TAPI E 10. ST7E 

TINE 1 

YDOT NFO 

Y ne o 

LAN 0 OA NLAN 

U NU 

HX, HY, H7, 

PX, PY, P7 6*N° 

TOTAL ANGULAR 
MOMENTUM VECTOR 
COMPONENTS <X,Y,Z). 3 

TOTAL l INEA F 
MC U £NTUM VECTOC 
COMPONENTS (X v V f '»). 3 

son y kinetic 

ENERGIES. NP 

RCOY POTENTIAL 
ENE D GIFS. NQ 


TOTAL A N CUL A R NC M ., 
TOTAL L INEAR HCM. , 


TOTAL K.E., 

TOTAL P.E., 

TOT AL ENERGY 5 

TIME 1 

YDOT NE 0 

Y NE 0 


TNis-RTIAL 
F* AMP . 


IV-41 



n n o o (i o .1 1 n n o n o o fl o o n n n .1 fl »i n .1 o o o o a n n d ci o o o n :i o M o .i o fl 


rfilMMT.cc^aT = 1615) ( JVPL ( J) , J = l, J C L > 


JVPLtJ) = INTEGER OENOTIHG GICPAL 
LOCATION CP JTH SELECTFO 
Va RTAELE FROM THE NTPLOT 
l C K 'G ARRAY, 

20 CONTINUE 


®EAO(MT, FORMAT = 575) NCI* ( NCO <T ) , T = 1 » 3) , MO-PID 

NCI = ELEMENT L CC ATI ON (LOCAL WPT 
JVPL ARRAY) PCP THF INCEFEN- 
DENT PLOT VARIABLE. 

NCC = ELEMENT LOCATION (LCOAL «pt 
JVFL AP°AY) FOR UP TO 1 
DEPENDENT variables TC plot 
^I f-ULTANECUSLV VERSES THE 
NCI OEPENTENT V A PT A p l. E . 

NGRin = NC. OP FLCT FRAMES TC USF 
F CP PLOTTING THF NO T - NCO 
GROUP. IE * THE NO. OF FRAMES 
TO USE SIOE °Y STOE TO 
EXHAUST THE RANGE CF THE 
INOEOENCENT VARIA°LE. 


IF (NCI .EC. t) GO TO 10CC 

THIS IS CUE TO PROCEOE TO NEXT SET 

THIS OPEN EN3EO LOOP PERMITS MANY SELECTION^ OF THE 
JVPL OATA VARIABLES WITH REGARD TO TNOF F ENOFNT ANn 
DEPENDENT VARIABLES IN ORDER TO FORK CROSS-PLOTS. 
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n n o .t o .i n 


t* 


c 

G 

G 

C 

r 

c 

G 

G 

C 

C 

c 

c 

c 

c 

G 

G 


-REAO(NIT,FCR*'AT = Afl, 2X, A8, ?X,6Afl> TITII, TTTID, 

* (PTITKII t I=l, 6 ) 

ALFMANUHFRin TITLING INFCRHA TICK 
TO INCLUDE ON PLOTTE 0 OUTPUT, 

TITLI = INDEPENDENT VARIABLE TlUr. 

TIUO = DEPENDENT VARIABLE TTTL r . 

PTITt = C VE° All TITLE ECP PARTICULAR 

identification CF this frame. 


GO TO 2C 
1C00 CONTINUE 

pc turn 

END 
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4 d001 N3dO OOfliSd <. = 
dOOl G3S013 9 = 
t iKU+Tj / i-li)) dOOl G3S010 S =■ 

(OH) aOOl 43d0 *r = 

(HO) dOOl M3dQ £ - 
(H) dSllQdiNOO Z - 
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* 5dAi N OliONOd ddiSNtfdi 

dHi ONilAdiSSmO d3031N i --- ddA.il = (P 4 1 )Adl 

- 310A0 Hir 3H1 dOi 

40Il«inwIS SI Hi dQd . 
dBOISNOO Oi SdlOAJ NOilUNDd 

d3dSNfldi BlVdtf d3S dJ * ON = 0AD4 

VltfQ lOoiNOO SiaAlsNtf 
AOSdHOddd ONlNlidU JA04 Ad & 3ZIS XldiVW 

(dX *6 ‘ UAUN * o * Adi ) HiQ«3d llvO- 


Nall3d ( H*» *ba* HV4l) d I 

NdfU3d (3HUH«J *0?* dl 

AOniS NIWWOG A0NdH03dd 

oazid^dMii aod 03 ddHn = 

3SN0dS3s 3HI1 

OiildVdNll dJd 3WIH^ - WVN1 

HtfNl (♦>» = iyAdOd^lNiGVjd 


*********************************************************************** 

* 
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c 

r 

c 

o 

r. 

c 

c 

c 

c 

r 

r 

0 

c 

c 

r 

r 

r 

0 

C 

C 

c 

r 

C 

c 

c 


r, 

c 

r 

C 

o 

c 

c 

c 

c 

n 

c 

c 

c 

c 

c 

r, 

r 

r 

C 

c 

r 

0 

C 

r 

c 

r 


2,j) = ITriN --- TR/VK5P£P FUNCTION r^ c t)T 

VARIABLE IDENTIFICATION. 
THIS T NT- GE P IS A LOCAL 
TOENTIFICATI ON IM;G C ®, 
REFERENCING (DE C E NOINT, 

U°ON ITYPE) PITHCR a 
SFNSOR SIGNAL CR A 
CCNTECLl'R CL' TP L' I VARI a"LE 
WHICH I? THE V (IN) OP 
THE EYPOrpSlCN -- 

V f OUT I 

---------- = TP 

V( TV) 


LRY{ 3 , J ) = JTFOUT 


transfer FUNCTION CUTOUT 
VARIABLE IDENTIFICATION*. 
THIS INTEGER I« A LOCAL 
I CE N T IFIC AT TON INTEGER, 
REFERENCING CCFCFkiOthG 
UFCN ITYPP» EITHER A 
SENIOR C TGN AL CP A 
CONTROL LE D CUTPljT VARIABLE 
WHICH IS TH c V l CUT ) CF 
THF EXPOESSTCN -- 


V10U T > 

VfIN> 


TP 


LRY( 4,J) = kOLOT 

KRLOT s j NO PLOTS 

KPLOT = 1 PLOTS WILL P- MACE 

LPY( E,J) =■ IAFLG 

this integer control parameter 

FERMITc THE USER TC SEl e CT 
THE CH ACACTERIRTT 0 BOOH FCR A 
GIVEN TRANSFER FUNCTION FRO 
ETTHFR THE OH AR ACT ER IS T TC H A TP TX , 
AR, OR ITS TP ANRCOSE , 

IAFLG = 1 PPCGPAN WILL USE CCTS 
FROM AP TPANSPCSF. 

OEFAULT VALUE IS 0 AND 
ROOTS FRCH AF WILL PE USED. 
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u o c? c 


0 

c 

c 

c 

c 

c 

r 

r 

r. 


r 

r. 

c 

c 

c 

c 

C 

c 

c 

c 

c 

r 

C' 

r 

C 

r 

c 

c 

r 

C 

C 

r 

c 

c 

c 

c 


NCTF - PROGRAM EXTRACTS RCOTS FOR ° OT H 

** AND ITS TRANSPOSE. THIS SERVES 
AS A SORT CF SFIF CHECK CN THE 
ROOT OIJALTTY. ALTHOUGH IT IS A 
RA D E ncruS°ANC c t SCOTS SRCK AP 
TRANSPOSE CAN* FE -CLEANER- than 
THOSE OFTATN' c O FROM AF. 

LRY( 6»J) = NC. Q V A RT A°L ES TC FEET '’AC* -- TTYRR 
H A X oc i r \/ARIA a LES CAN P£ FE n ° AC K FOR 
THE TYRE 7 RSELCC OPEN LCCR 
TRANSFER FUNCTION. 


L R Y ( 

7,J) = LOCAL 

ID. 

Cf 

FIRST 

n 

TC 

RETAIN 

LRY ( 

*,J) = LOCAL 

in. 

OF 

SECONO 

R 

TC 

RETAIN 

L R Y { 

<3,J> = LOCAL 

TC. 

CF 

THIRO 

0 

TC 

retain 

CALL RfacIk 

(IPY, 3, NCYC 

* -3, 

KR) 





MATRIX SITE 3 °V NCYC OEFTM NG EXPONENT PC c 
TOLERANCES TOL = <10.»*«EXP 

FOR THE JTH CYCLE - 

IRY(l.J) = ROOT TOLERANCE EXPONENT 
IRY(2,J) = GAIN TOLFOanCE EXPONENT 
IRV(3,J) •= ccot tolerance exponent USEr TC 

REHOVE SHIFT FREQUENCY (SU°POUTTNP NIJHS) 


NOTE - I? ROOT OR GAIN LE TOL, SET 
ROOT OR GAIN EQ 0. 
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C 

C 

0 

c 

c 

c 

r. 

c 

r. 

c 

c 

c 

c 

c 

c 

c 

r 

” 9 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c. 

0 

c 

c 

c 


no 500 icvc=i,ncyc 

IF ( I TYPE ,E0. 5) GC TO 5G0 

•PFflDf NIT t PCPYflT = 20A4) ( TITLFC I) ,1^1,20) 

flO C HAPACTER TITLE FOR TRANSFER FUNCTION 
ITENTIFir ATICN 

FEA0(NIT, FORMAT = 5A4) <L PNAME (II , 1 = 1 » 5> 

IFNAMEU) PERMITS UP TO 5 FOUR CHARACTER 

IOENT IF ICATTONS WHICH SELECT THE 
PICT niSFlAY MCCE. 

LRKAMECT) = 4 H {AIL °L ANK) NO 0XS c lA*S 

ARE IMPLE m ENTEO GO TO 500 

THE CH ACAOTEF 1ST IC PCCTc 
C C° The SYSTEM ARE FCUNP. 

LF N A*E III = 4 HP OP r ONLY A POPE DISPLAY, 

= 4HNT0V ONLY A NICHOLS DISPLAY. 

= 4HNYOU CM Y A NVOUTST DISPLAY, 

= 4WN T MV °CTH NICHOLS A NO NY0U T ST, 

= 4HPCMN GIVES p O n E» NTOHOLS, NYOUTST, 

- 4 HROO T GIVES A ROOT LOCUS OT^RIAY, 
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o r> n n n o n ,*1 n r> o o r> n r> r> ri o n n ^ n n no i o o o o o n o o o o o o 


r 


no 5G0 I C F = 1 » 5 


TF (L p N A M F (TOP I .EP. 4H > GC TO FOC 

IF (LPNflMF(ICO) ,EO. 4HDC0F 

*.OP. L°Kfl^EtICP> • EO . 4HMCH 

*.0R. LPKAVPUCCI ,EQ. 4HK Y OU 

•.C^. L PK AHt (TOPI .cQ. ^MNV 

*.O c . LPKAME(ICF) .EO. 4HFCNM GC TO ?0S 

IF (LPNA^E(IOP) .£0, AHFOOT) GO TO 300 


2*0 OCNTINUP 


»*» FPEOUENCY FESPCNSE SECTION *** 


FE801MT, FCP^AT = 6 F 1 C . G > F>!n, F*AX, OBMTK, 0P>1AV* A m T*, «“fiv 

F w I K r FFPfJUEK'CY FW£CF L CWEP 1 1 KT T 
f^av - FDFQLFNrV FWFpp UP C FR LIMT 

OPvjh = "TMKUM O’ AHFLTTICE for BOCE, NICHOLS *>1 CTS 

ORHAX = ‘'AYT^JM or amplituof for "ooe» MCHCLS PLOTS 

A*IN = KINIKJM AMFLITLiOE FOP N YOU 1ST PLOT* 

A H A X = *AXlHUH AMFLITUOE FOR hYCUIST FL CTS 


GO TO FOG 


3G0 CONTINUE 
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c 

r 

c 

c 

r 

C 

c 

c 

c 

0 

C 

r 

c 

c 


r 

r* 

i 

c 

c 

o 

c 

/ 

o 

r 

c 

c 

r 

r 

r 

C 

r 

C 

c 

c 

r 

c 

r 

c 

c 

c 

c 

r 

c 

£ 

c 

r 

c 

t 

c 

c 


M * RCCT ICCUS SECTTCN *** 


CALL C EA0IM (IJM, 2, NRLC» ?, KRJ 

MATRIX SI 7E 2 °Y VRLC FOR RCCT LCCHS PICT CONTROL 

NRLT r NUYRCR CF PCCT LOCI TC PFPFCRV 

C CP THE JTH ROOT l CCT - 

= ISNIM = l STARTING PC I NT T c 
OPEN L CCP ZFPC. 

= Z STAPLING PCTN'T 79 
0 °t N l CCP FCLE. 

= t starting pctnt re 

CLOSED L CCP P CL F • 


S ELEMENT LOCATION in OfCT ARP AY 
F C° STARTING PCOT l CCT, 


CALL PE AO (m, e, NRLC, YR, YR) 

MATRIX SIZE 6 BY MfiLC FOP RCCT ICCUS CONTROL CATA 
FOP TH^ JTH RCCT L CCT - 

W1(1,J> s T MET A C ( J) INITIAL SF ARCH ANC-tr 

NCR MA ll Y -1«C. (DEGREES) 

W1(2,J) = eCL SCALE FACTOR, NC^ALIY 
SCI = l.G 

Wi(’,J> ALOC PHASE CCNTRPL PARAMETER. 

AlCC = *1, --- J.AC, OEG. PHASE. 
ALOC = -1, 0 CRC. PHASE. 

= Y MI N MTN REAL VALUE TC PLOT. 

W1(E,J) = XMAX MAX REAL VALUE T 0 PLOT, 

HKf,J» = YMAX MAX I HAG VALUE TO FLCT AMO 

YMTN SET TO - Yvax. 


5 C 3 CCNTINUE 

RETU°N 
r NO 
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V. 


PROGRAM OUTPUTS 


This section discusses the various program output information 
and correlates the output data with both the input data and the 
problem simulation. This information is presented in much the 
same fashion as was the input data stream of the previous chap- 
ter so as to better acquaint the reader with the actual formated 
output as it is presented by the program. 

It is pointed out that the output stream will not reflect cer- 
tain outputs that occur from routines that identify troublesome 
areas such as matrix singularities. Recall also that the basic 
input routines READ and READIM can also print out input matrix 
data as dictated by the user. These printouts will not be in- 
cluded either. Reference is made to the theoretical volume 
(Vol. I) and to the input data stream (Chapter IV) to correlate 
certain outputs with both the theory and the user input require- 
ments . 
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THE FOLLOWING L 7 ST TnFNTlFIFS TMf 
OUTPUT VARIABLES SUMMARIZED ON ThF 
PREVIOUS PAGE. 


C 

c 
r 
c 
c 
c 
c 
c 
c 
c 

r 

c 


r 

c 

NR = 

MO. OF 

POOJES. 



c 

f 

NH = 

MO. OF 

HINGES. 



c 

c 

NSPT = 

MO. OF 

SFNSOP POINTS* 



c 

f 

NOFLTA = 

MO. OF 

CONTROL SYSTEM OF| TAP. 

c 

r 

NU r 

MO. OF 

U*S. 



c 

o 

NBFTa = 

MO. OF 

RFTa*P. 



0 

r 

NL AM s 

MO. OF 

LAMBDA**? (CONSTRa JNTP) . 

r 

c 

MFO s 

MO. OF 

STATF fOUAT I O mS . 



c 

c 

NPMAX = 

MAXIMUM 

OIMFMSIONFO NO. 

OF 

ROOTFP. 

c 

o 

NHMAX = 

MAXIMUM 

oimfnstonfo NO. 

OF 

W T NOES * 

c 

c 

nPPmax = 

MAXIMUM 

oimfnsionfd no. 

OF 

SFNSOP POTNTS. 

c 

0 

MMWMAX s 

MAXIMUM 

OIMFNSIONEO no. 

OF 

MOM. WHFFj S. 

c 

MMWROD = 

MAXIMUM 

niMFNSIONED NO. 

OF 

MOM. WHFFj $ 

0 

c 





PFR ROOY. 

c 

c 

nmoboo = 

MAXIMUM 

DIMENSIONED NO. 

OF 

MOOFP PFR BODY 

c 

KMU s 

MAXIMUM 

OIMFNSIONFO no. 

OF 

U*S PFR ROOY, 

r 

r 

s 

* ♦ NMOROO ♦ NMWBOO. 



V 

c 

c 

KV r 

MAXIMUM 

OTMFNSIONEO PT7F 

FOP STATF vfctop 

c 

KU r 

MAXIMUM 

OIMFNSIONEO NO. 

OF 

u»s* 

c 

= 

MRMAX*(* ♦ NMOBOO) * NMWMAX 

• 


c 

c 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

r 

c 

c 

c 

c 

c 

r 

c 

c 

r 

c 

c 

c 

c 

c 

c 

c 

c 

r 

c 

c 

c 

r 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


OUTPUT VARIABLE IOFNTTFIC ATION SUMMARY 
(CONTtn) 




STARTT = START TTMF FOP TIME PESPONSF. 
OFLTaT SS integration stfp size. 

FNDT a FK|r> TIMF PGR TIME RESRONSF. 


01 = X COMPONENT OF GRAVITY VFCTOR. (INPUT) 

G2 c Y COMPONENT OF GRAVITY VFCTO&. (INPUT) 

G3 s 7 COMPONENT OF GRAVITY VECTOR. (INPUT) 

GMAG = SQPT(G1**2 ♦ G2**2 ♦ G3**2) — ACC. OF GRAVITY. 

GAM1 = OIPECTTON COStNE (GRAVITY VFCTOR ANO X). 

GAM2 « OIPECTTON COSINE (GRAVITY VFCTOR ANO Y) . 

GAM3 = OIPECTTON COSINE (GRAVITY VFCTOR ANO 7). 

PCMAG « PFFEPENCE RADIUS FOR ACTING GRAVITY VFCTOR. (INPUT) 

NOPRNT s MULTIPLE OF nELTAT TO PRINT time RESPONSE. 

NOPLOT » MULTIPLE OF OFLTAT TO WRITE PICT TAPE. 

TFLNER s 0 (NON LINFAR TIME RESPONSE) 

a 1 (LINEAR ANALYSIS OR FPFO RESPONSE) 
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c 

c 

c 

r 

c 

C THE: FOLLOWING SF&IFS OF ARRAYS PFPT&TM 

C TO THF MODEL SIMULATION! AND ARE PR TNTFO 

r OUT AT Thf REGINNING OF EACH SIMULATION 

c 


c 

r 

C - — THF TOPOLOGY ARRAY (ITOPOD FOR THIS CASE FOLI OWS 

f 

C (1) (?)... (WH) 

r 

C — - 1 1 

c 2 1 

c 

C THIS IS THF INPUT I NTFGFR APPAY TTOPOl 

C (RFF PREVIOUS CHAPTER — INPUT OaTA) 

c 

c 

c THF CONSTRAINT SorrtFTCATTONS FOR THIS CASF FOU OW 

C 

c ( 1 ) (9) ... (NH) 

C 

C — l 1 

C 2 ) 

o — — — — «— 3 l This is thf input intfgfr appay thoata 

r — 4 1 (SFF PREVIOUS CHAPTER — INPUT OATA) 

c S 1 

c 6 1 

c 7 1 

c 

c 

c — - — Tw£ SPECIFIED INITIAL HINGE ANGLFS 

C AMO DISPLACEMENTS (RETAH) FOLLOW 

C 

c U) {?) ... (NH) 

c 

C - — 1 1 

c — -- — -- 2 1 ROWS 1-3 = HINGF ANGLES (CONST STFNT WITH TTYPE). 

r — 

C- . . ROWR 4 -6 = HINGE 0 I SPL A CEMENTS- 0 RFf ATTVF TO P # 

C 

c G 1 

C 

C 


« 
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C 

c 

c 

r 

c 

c TWF SorclFTFO INTTaL wINGF RATFS (BETAHO) FOLIOW 

c 

f (1) (?) ... (NH) 

c 

C — -- 1 1 

C — 2 1 ROWS 1-3 = ANGULAR RATFS. 

C — 

C-< — — . . ROWS 4-6 = DISPLACEMENT PATFS- 0 RF| attvf to p. 

C— ———————— . . 

C — A 1 

c 

<■ . . . i 

c-- the no. or flaSTtc moreS/rody array (IRgflxj foi lows j 

C ! 

C (1> (?) ... (NR) I 

C - ; : . ... | 

C — — 1 1 THE JTH FNTRY IS THF NO. OF FLASTTC I 

C MOOES FOR OOOY j. 

C 

r 

C— — THF NO. of p/ 0 HTNOF oOINTS/ROOY ARRAY (NHPOT) FOLLOWS 

r i 

C (1) (?) ... (NR) 

r 

c — — 1 1 Ft FMFNTS are ThF NO. of p/0 

C HTNGF POINTS On EACH PODY, 

c I 

c I 

c TMF NO. OF SENSOR POlwTS/ROOY ARRAY (NSPOT) FOl I OWS • 

<* i 

C (1) (?) ... (NR) 1 

C 

r 1 1 Ft FMFNTS are The no. of sensor 

C POTNTS ON EACH PODY. • 

C 

C 

C ——THF MOM. WHEFL/ROOY TABLE ( N M OW j FOLLOWS 

C 

C (1) (?) ... (NR) 

C 

C — 1 1 COL J = JTW BODY 

c — — — 2 1 ROW 1 * NO. OF mom, WHEELS ON RORY J. 


. . ROW ? = NO. OF VARIARLF SPEFO WhFFI.S ON ROOY J. 

- . ROW 3 = SUCCESSIVE ROWS ARF THE MOMENTUM WHEEL NUMBERS 

- . . . ON RODY J (IN ASCENDING OROFP). 

p.nmwroo 1 . 

c 

c 

c 
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THF STaTE VECTOR L FNGTH ARRAY (lENU) FOLLOW*; 

(1) (?) ... <?*** ♦ ?> 

1 1 

THE ELEMENTS ape the i.fngths OF segments 

OF THE FTATE VECTOR. 

ORDER TS 

‘. 1(1 ) *U( 2 > ♦,.»U(NR) »XT ( 1 ) •..♦XT (NR) .RET A. OF) TA 

THE STaTE VECTOR LOCATION ARRAY (LOCU) FOLLOWS 
(1) (?) ... (?*NR ♦ 2 ) 


1 1 

LEADING FLFMENT LOCATION in THE STATE VFCTOR 
FOR The SEGMENTS DESCRIBED TM ARRAY (lENU). 


THE SPECIEIED SENSOR POTNT/ROOY CORRELATION ARRAY (TFTSmW) EOLLOWC 

(1) (?) (NSPT) 

1 THE JTH ELE m FNT TS THE BODY NO. ON 
WHICH SENSOR POTNT v J IS LOCATED* 


°^JY?! 4 Zp 
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r 

c 

r 

c 

c 

r TMF FOLLOWING data TS SPFCTFTFn MOM. WHFEL INFORMATION (IF ANY) 

f AND CONTROL L p P INFORMATION 

C : 

C 

c - THF SPFCIFIFO mom. WHFFL CONTROL array (IMO) roi I.OWC 

r 

r (1) (?) ... ( NOFMO ) 

c 

r ■ COL .1 = JTH MOMPNTUM whffl 

1 1 ROW 1 = WHFFL SPNSOR point no. 

2 1 ROW ? = SPIN AXTS 

3 I ROW 3=1 ACTIVF 

= 0 CONSTANT speed 

r 

o -»TWF SPFCIFIFO MOM. WHFFL RaTFS aNO INFRTIAS (AMO) FOLI OW 

C 

c (1) (?) ... (NOFMO) 

C 

C COL J = JTH MOMFNTUM WHFFL 

C 1 1 ROW 1 = INITIAL WHFFL SPIN RATF 

r ? 1 ROW ? = SPIN IMFRTTA 

r 

C ’ 

r-- —The SPFCIFIEO CONTROLLER TNITTAL conditions 

c ANO CHARACTERISTICS FOLLOW 

c (THF FIRST NDFl T a A°F INITIAL CONTROLLER STaTF 

r VARIABLES* THFRF are K AnnTTIONAL PARAMETERS) 

C 

C THIS IS THF USER INPUT ARRAY CNTOTA. 

C 

C THF ADOTTTOMAL K PARAMETERS » TF ANY. ARF 

C AVAIL ARLF TO THF USER FOR USER-PAK DATA. 

C ... 

C THE FIRST NOELT A ENTRIES IN THIS ARRAY 

C ARE THE INITIAL CONDITIONS FOP THF 

C CONTROL VARTABLFS. 

C 

C 

C 

r ********************************************************************** 



* 
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noon 


c********************" 0 ** ** ** **** * 


C* GUBPOUTINF MRIfilO — OUTPUT FOR RTGTD POOY 


******* ******************** ** 

THF following IS TYPICAL for thf ith rooy 


C(jmm ARY OF input OATa fop POOY T WHICH T< RTGIO. 


c .-thf 6 X 6 infrtta matrix is — 

c 


1 

? 

1 

A 

G 

6 


( 1 ) 

1 I XX 
1 -TYX 
1 -T7X 
0 

1 -S 7 

1 SY 


1 


( 2 ) 

• I X Y 
TYY 
-IZY 
S 7 
0 

-sx 


( 3 ) 

- T X 7 
-TY 7 
TZ 7 
-GY 
SX 
0 


( 4 ) 

0 

S7 

-SY 

M 

n 

n 


<*> 

-sz 

0 

sx 

0 

M 

0 


( 6 ) 

SY 

-sx 

0 

n 

o 

M 


r FOR POOY I THE P -0 HINGE NO. ANn THE EULER ROTATION TYQF 

r ftppFAP tn the Foil owing Integer array which ts foli owfo pv 

K w ^ ^ 1 _ * _ « . V .km nArTTTniu 


c 

r 

c 

c 

c 

c 

c 

c 

c 

c 

c 


AN APPAY CONTAINING FULFP ANGLFS ANn 

vfctop COMPONFNTS ( 4 ♦ S * 6 ) THAT POSITION Thf HTMGF TPIAO 

WPT THE POOY TRIAD 

DATA HERE 


TF POOY T HAG ANY SFNSOR POINTS. 
THE FOLLOWING WILL PE PRINTED 


___FOR BODY T the SENSOR POTNT NO. AND THF EIJIFR POTATION TYPE 

r APPEAR IN THE FOl LOWING INTEGER ARRAY WHICH TS FOILOWFO 


c 

C 

c 

c 

c 

c 


BY AN APPAY CONTAINING EUl.FR ANGLES AND POSITION 

VFCTOP COMPONENTS ( 4 .S.G) THAT POSITION THE SFNSOR TRIAD 
wrt the POOY TRIAD 

DATA HERE — — 


• .UK******************************************************************!! 
********************************************************************* 


V-9 



c 

C»*«-»«*********** a********************** ******************************** 

c* 

c* qnoooilTINF MSMOnC - OUTPUT FOP FLFXTRI.F BODY, CONSTSTFNT MA«S MATRIX 

r* 

C*********************************************************************** 

c 

c 

c TYPICAL OUTPUT FOP TTH BODY 

r 

C 

c THIS SUBROUTINE PRINTS out several matrices THAT apf 

c RFLATFO to THF FORM of thf governing fouattons. 

C - , . .. . 

r the reader will of referred to vol T and thf tnput 

c OATA FTRFAM for fijrthfr CLARIFICATION. 

c 

c ; SIJMMAPY of INPUT oata fob rooy t which ts 

C Fj EXIBl.F w/consiftfnt mass matptx. 

c 

c 

c— —thf input parameters Ifpbm, tfoiak. tfdiad apf 

c 

c RrE TNPUT DATA FOR MSMODC 

c 

c 

c— THE JOOF TaBLF FOU OWF— *- 

c 

c nEGPFE OF FRFFO0M as tnput 

c 

C FFF TNPUT DATA FOR MSMOOC 

c 

c 

c — —THF MODE SELECTION VECTOR FOLLOWS 

c 

c MOOF SELECTION VFCTOR AS TNPUT 

C SFF INPUT DATA FOR mSmoDC 

c 

C -- -- - • 

c .—FOP ROOY NO. I THF POSITION VFCTOP FROM THF ROnY ORTGTN 

c TO JOINT K IS 

c 

C X = Y = 7 s 

C 

c WMFRF K IS JOINT COOPDINATFS 

c USED TO OFVELOP RIGID RODY MODFS 

C 

r FFF INPUT DATA FOR MSMODC 

C 



n n n n n m ^ ° 



C 

C 

C 

C 

r 

C- 

C 

C 

C 

r 

c 

c 

c- 

r 

c 

c 

c 

c 

r 


TWF COMSTSTFNTt DFPftRTITinMFn maSS MaTRTX TS 

TMT^ is THF PFOAPTITIOMFO 
matrtX ANn T s rOMSTSTFNT WITH 
ThF IDOF TARLF. 


>Th f RFPARTTTIONFn won/'L MATRIX TS 


THTR is THF RFPARTITTOMFn Monai 

matrix, thf °ows arf conststfnt 
with THE RFPAPTTTinNFP M *S«; matrix 
A ign thf POLS ARF CONSIFTFMT WTTH 
THF FLE m FnTF OF THF MOPF SFI Pf T TON 
vr r T op . 


TW F -UN n FFORMFO~ TNFPTTA MATOTX (MU) TS 

THTS IS thf mo matrix mOTFH as 
EQUATION 11-87 (VOL I) • 


f- 

r 

c 

r 

c 

c 

r 

r 

r 

c 

r 

c 

r 

r 

c- 

r 

r 

r 

r 

c 

r 

c 

r 


jMPPfr thfn fou nwq MATRTres 

A COFFF T r IFMT^ 

q COFFFiriFMTS 

rorxv COrFF ICIFNTS 
cnr X7 rOFFFirTFMTS 

rnF Y7 rOFFFir tfnts 

WHICH &PF thf Al’PMAf Rf AND r COFFF I C I FNJTS 
GTVFKJ in THF F X PRF^c; I ON 

EQUATION 11-88 (VOL I) 


JMFQF THFNJ FOLl.HWS maTRTCF^ 

rn * C??* CT?, 

n* f ci3t c?3 

W HirH ARF inFNTTFTFn I M 

EQUATION 11-89 (VOL I) 


* 



V-ll 



n-A 


« 


»««*»« 


— — da3H VJLVU - 

UVIdi AUUb iHl ldrt (JVldi dOSNiS 
dHl NOillSOa 1VH1 SJHONV di 1H3 UNlNiVlNOO Atfdd* NV Ad 
UiAOllOi SI H?IHM AVddV d3t>31NI 9Nl«0ll0d 3Hi Nl dV3ddV 
iNXOd dUSN3S 3Hi 01 9NKlN0dS3ddUU *0N IN I OP iHl UNV 
ddAi NOilVIUd d3in3 3H1 4 # ON iNIOd dOSNdS 3H1 I AUOb dOJ 


CiilNldd id Til* 9NI«0n."»0i iHl 
4 SiMOd dQbNiS ANV SVH I AUOd di 


— — idiH VlVU - — - 

UVldl AUUd dHl Ad* GVldl 
idNlrt iHl NOl liSOd 1VN1 Si^UNV di lOd 9NiNi*ANOO A^ddV 
NV Ab UirtOl lOd HOIrtW AVddV db93lNl 9NIWU1 IOi iHl Nl dVdadV 
“ * 39N1H 0~d 3Hi 01 UNlUNOdSiddUJ *0N INlQP iHl ONV 

idAl NOilVlOd dilOi JHi ‘•UN 39NlH~0-d 3H1 I AUUb dOi 


•AliAllUidSid SiIllOO"liA IVQON “IViliM 
UNV SANiNiOVHdSlU iVUOw 1V111NI 

Jh 1 9N I N 1 V 1N00 SAVddV 0*1 S*Q~HUi idiHi 


idiH VIVO 

— bl XldlVh UNIdWVO 1VUUN iHl 

itiiH V1VU 

— bl bbiNJiUS 1VU0W iHl 


**U 

*♦0 

0 

J 

0 

3 

3 

3 

J 

3 

3 

— 3 
3 
3 
3 
3 
3 
J 
3 
J 
3 
3 
3 
J 

■ — 3 

3 

3 

J 

3 

3 

3 

• — 3 
3 

J 

3 

J 

— -J 

3 

3 

3 

3 

- — 3 

3 

3 

3 

3 

3 

J 



****************** 


<hhhhhmmhmmmhkhhh>*** # * 


r 
r 
c 
c 
r 

c a******************************"''”" 

c' .iioqoiiTTNF M5MO0L - OUTPUT FOP FLF«,RLF ROOY. <-U»°FP m 4 FS «iT»TX 

"... * 

r 
r 
c 
r 
r 
c 
r 
r 
r 
c 
c 
c 

£ ..output matrices tnfro. stato, masso, oocoff, aocoff, fpoff 

r ARC THF I. -S. M, 0 . A, F. P 4 PTTTTONF 

j; . Sfspfctivfly of the matrix mo of equation 11-87 

r TN VOL T . 

r 


TYPICAL OUTPUT FO B TmE ITH RODY 


TwT^ SUPPOUTTNF PPTMTP OUT SFVFPAL MATRICES THAT APF 
PFLATFO TO THF form of thf oovfrninr FOLIATIONS. 

TMF PFaOFP WILL RE PFFERRFn TO VO| I AMO THF INPUT 
n*TA CTRFAM FOR furthfP Cl ARIET oaTION. 


C 

c- 

r 

c 

c- 

r 

C 

c 

c 

c 

c 


OUTPUT 


OUTPUT 


MATRIX 

matrices 


MUO TS haTRTX mo of EQUATION II -87 (VOL I) . 


A COF . RCOF* CXY, CXZ* CYZ« 

APF THF alpha. R. ANO c 
roFFFlCIFMTS TN THE maTPIX 
OTVFM AS EQUATION 11-88 (VOL I) 


« 
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c 

c 

c 

c 

c 

c 

c 

c —OUTPUT MATRICES CUt C22. C33. 

r c i?. c 1 3 « c? 1 

c 

c ARF THE CONSTITUENTS OF THE 

C MaTRTX STVEN AS EQUATION 11-89 (VOL I) 

C 

c 

c —OUTPUT MATRIX XFO 

C 

c contains thf initial modal 

c OFFLFCTTONS. (AS INPUT) 

c 

c 

o — — — OUTPUT MATRIX XEnn 
C 

C contains thf INITIAL modal 

c VFlOrlTlFS. (AS INPUT) 

c 

c 

0— — — FOR RODY I THE P-0 hImGF no.. ThF FULFR ROTaTTON TYPE 
C AND Thf JOTNT NO. CORRESPONDING TO THF P-0 HIMGF 

c APPEAR IN THE FOLLOWING INTFGFR ARRAY WHICH is FOLlOWFn BY AN 

C ARRAY CONTAINING Fill FR ANGLFS ThaT POSITION THF HTNGF 

C TRIAD WRT THE BODY TRTAO 

C . _ 

C DATA HERE 

C * - . . , .... , 

c 

c IF RODY I HAS AMY SFNSOR POINTS, 

c thf following will rf printed 

c 
c 

c _; -__eor body i thf SFNSOR point no., the euifr rotation typf 

c and Thf joint NO. CORRESPONDING to thf sensor potnt 

c APPFAR IN THE FOLLOWING INTEGFR ARRAY WHICH IS FOLI OWFD 

C PY AN ARRAY CONTAINING FULFR ANGLFS THAT POSITION THE 

C SFNSOR TRIAO WRT THF RODY TRIAD 

C . . 

C — - oaTa herf 

c 

c 

c 

p**********#***********************************************"*********** 
* 


V-14 



c 

c 

c 

c 

c 

c 

c 

c 

r 

c 

c« 

c 

c 

c 

c 

c 

c 

c 

c 

r 

c 

c 

r 

c 

c 

c 

c 

r 

c 

c 

c 

c 

c 

r 

c 

c 

c 

r 


IF THF USER UTILT7ES POLYNOMIAL INPUT FOR 
CONTROL SYSTEM TRANSFER FUNCTIONS* the following 
MATRIX WILL RE PRINTED, 


OUTPUT MATRIX CRLY <*P| Y*?#NPLY) 

where KPL Y = ROW DIMENSION SIZE OF CPLY IN SUBROUTINE CONTR[ 

nply = no. of input polynomial Ratios 

ANn FOR I = ODD INTEGER — 

CO L I = DENOMINATOR POLYNOMIAL coefficients 
TM ASCENDING ORDER 

COL J*1 = NUMFRATOR POLYNOMIAL COFFFTCTFNTS 
IN ASCENDING ORDER 


THE F0| LOWING INTFGFR ARRAY (TNnEP) PRFSCRTRFS INDEPENDENT 
VARTAPI.FR ( 1 ) AND DEPENDENT VAPIARI.ES (0) 

THF FLEMFNTS OF THIS ARRAY IDENTIFY 
WH I CM V AP I AR|_ F S SURVIVE IN THF FINDU 
SEARCH To dftfrminf AN INDFPFNDFNT 
SET TO °F INTEGRATED. 

Tme SURVIVING VARTARtES WILL ALSO 
pppRrSFNT THE FIRST NX ROWS OF 
THF LINEARIZED MATRIX* A, USFn 
PY SURPOUTTNF DYNS40. 

MOTE - NX = NO. OF NON-ZFRO ENTRTTS TM ARRAY TMoEP 


C 

C 


C 

c 

c 


* 
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c 

c 

c 

r 

c 

£#»«««*««»»«•«*#«#**«#+***«««*»**+***************************«********** 

c* 

c* c^noonuTINF OYNS?0 OUTPUT 
C* 

C* #**»*#* ********************************************************* 

r 

c 

C THF PPINTOUT TS TYPICAL FOP A GIVEN SIMULATION TTMF, T. 

r 

r THF TrO PRINT OUT IS ALWAYS GIVEN 

c (FVFKI TOR A { TNFAPT7ED ANALYST ^ ) 

r 

r 

c the data arf ore SENTfo IN VECTOR FORM ano 

c ORnFRFD as roi Vows — - 

c 

C — THF STATE VECTOR V s 

c 

c -.-THF STATE VECTOR TTmF DERIVATIVE YDT * 

c 

c —THE BFTaS (EULFR ANOLFSf POSITION COOROTNATFS) ARF 

r 

C--> --THE RET A TIME DERIVATIVES ARE 

r 

c ......THF OFlTAS (C0NTOOI SYSTEM VAQTaRLES) ape 

c 

C-— — — THE OFl T A TIME HEP T VAT I VES ARE 

r 

c 

c 

* 
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THF: FOLLOWING INFORMATION is typical for rooy t 


c 

c 

c 

c 

c 

c 

c 

r 


c 

c— 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r 

c 

c 

c 

r 

c 

c 

r 

r 

c 

c 

r 

r 

c- 

c 

r 

c 

c 

c 

c 

c 

r 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


for ROoY I 

ORPFR IS 


FOP ROOY I 
OOOEP 


TMF VEi OCTTTES APF 

omox 

OMOY ROOY ANGULAR VELOCITY 

OMG7 

V ROOY TRANSLATIONAL VELOCITY 

w 

x t n ( n 

XTO (?) 

rooy moral velocities 


XTO (NF) 

THFTAOM) 

t MOMFNTUM WHEEL ANGULAR VFLOCTTY 

| (RELATIVE TO SENSOR POTNT TRTAO) 

(ACTIVE. WHEEL) 


the corresponding momenta are 

O0DY AXFS REF. ANGULAR MOMFNTUM 
(INCLS CONTRIBUTION OF CONST. SPEED WHEEL) 

ROOY AXES REF. LINEAR MOMENTUM 

) 

P XT (?) 

ROOY AXES PEE. MODAL MOMFNTUM 

P XT (NF) 

H mw ( 1 ) 

ROOY AXES REF. MOM WHEEL MOMFNTUM 

| (ACTTVE whfel> 


IS HX 
HY 
H7 

l x 

I.Y 

L7 
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o o o r> n -> 


C 

c 

c 

c 

c 

c 

c- 

c 

c 

r 

c 

c 

r 

c 

c 

c 

r 

c 

c- 

c 

c 

c 

c« 

c 

c 

c- 

r 

r 

c« 

c 

c 

c 

c 

c- 

c 

c 

c 

c 

c 

c 

c 


--FOR ROOY I ITS CONTRIBUTION Tn TOTAL ANGULAR 
A NO LTNFAP MOMrNTUM TS 

opoer is MX 

HY ANGULAR MOMENTUM 
M 7 

(RFFFRFNCFD to INERTIAL ORIGIN) 
LX 

LY LINEAR MOMENTUM 
L 7 

—FOR ROOY I ITS CONTRIBUTION TO TOTAL KINETIC ANO 
POTENTIAI energies IS 

— FOP ROOY I THF ELASTIC OFFLFCTTONS ARE 


THE INTERCONNECTION CONSTpaINT FORCES (I AMRO A^ ) ARE 


-tme total angular momentum vector is 

X f Y ♦ Z COMPONENTS in tnfpttal REE. axis system. 


TWF TOTAL LTNFAP MOMENTUM VFCTOP IS 

XtYfZ COMPONENTS IN INEPTI AL PEE. AXTS SYSTEM. 


the total angular MOMENTUM s 
twf Total linear momentum = 
THE TOTAL KINETIC ENERGY « 
Twf TOTAL POTENTIAL ENERGY = 
THE TOTAL ENERGY (T*V) * 


t #•«««««»« 


********************************** 


«« 

«« 


*** 


* 
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o o n 


C 

c •*»#***** *#***«# 

c* 

C* <;iiqoouTTNF 


0 VMc; 3 0 OUTPUT 


r* 

c 

r 




C 

c« 

c 

c 

r 

c 

c 

c 

c 

c 

r 

c 

r 

c 

c 

c 

c 


SUMMARY nr PLOTTING information 

THF OUTPUT SUMMARIZES THF input data 
which comtpollfp thf time response 
PLOT logic* 

SEE INPUT DATA FOP nVNS^O FOR 
OFT A T t FO nFSCP TPT T OnS . 


NOTF- THF OUTPUT VAPIARLF NPPLOT 

IS THF NO* OF LOGICAL PFCOPOS 
WPITTFM ON THF PLOT TAPP. 


*********************** 

*********************** 


******************************* ****** * 
*********************************************** 



r> d o o 


C 

c 

c 

c 

c 


c 

r * «»««««« **«♦♦******# ********************* 




c* 

C * cuPonuTTNF 0YNS4n OUTPUT - I INFapTZFD SYSTFM analyst?: 

c* 


c 

c 

c 

c 

c 

r 

r 

c 

c 

r 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r 

c 

c 

c 

c 

c 

c 

c 


OUTPUT MATRIX -A- (N.IOtNX) 

M JO = total no. OF EQUATIONS 
LTUFART7F0 (INCLUDING 
AUXILIARY EQUATIONS FPOM 

surpouttnf eq A nn > 

NX = NO OF TNOFPFNDFNT ST A TF EQUATIONS 
nETFPMTNFD RY FINOU IN OVNSZft. 


TWTS IS THF MATRIX OF PARTIAL DERIVATIVES 
WWTCM ARE THE LINEARIZED COMPONENTS OF THF 
STATE V AP i ARLES AS OETF RM INFO RY 
SURROUTINE ltneap. 


the OROEP OF THE VARTARLFS for THF 
ROWS IS 

SIZE 

PLANT VAPTARLFS NY 

control variablfs ndelta 
contpollfr outputs nbto 

SENSOR SIGNALS NXSS 


THF COLUMN ORDFR IS 

PLANT VARIABLES NY 

control variarles noelta 
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OUTPUT MATRIX -t- <nx.nx) 


c. 

c 

c 

c 

c 

c 

c- 

r 

c 

c 

c 

c 

c 

c 

r 

r 


ThTS IS IMF SIMILARITY 
transformation matrix 
THAT INTRODUCES THF (NbTO ♦ NXSS) 
AHXJl IAPY VARTaRLES INTO THF 

transformed state eouations. 

RTP. MATRIX R TN EQUATION III- 24 (VOL I) 


C 

c OUTPUT 

C 

C 

c 

c 

r 

c 

c 

C-- OUTPUT 

C 

c 

r 

c 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


matrix Y* <1*nX) 

this is Thf transformed sTatf 

VECTOR INITIAL CONDITIONS. 

RFF. VrCTOP 7 TN EQUATION III-24 (VOL I) 


MATRIX a* (NX, NX) 

This IS thf TRANSFORMED, LINFART7FD 
STATF VAPlARLE COEFFICIENT MATRIX 

THAT ts the basis pop thf entire 
linearization PACKAGE. 

pr E. EQUATION III-28 (VOL I) 

THF ROW/COL VARTAPLE ORDERING ANO 
SIZES ARE S 

VARIABLE TD. qT7r 


PI ANT VARIAniF*: 

NY? 

plant sensor stgnals 

NXSS 

CONTROL SYSTFM vartarlfs 

NO? 

CONTROL OUTPUT*; <R*S) 

N«TO 


NOTE- NY? s NY - NX SS 


NOR x NDFL T A - NBTQ 


o 
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PT A 


PTA* 


r 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


RF A| P APT IMAGINARY PART RFAL PART IMAGINARY PAPT 

COMP| FX ROOTS ORTA INFO FROM 
A AND A* PFSPECTIVELY. 

ThFSF aRF ThF POLFS of 
thf closed loop system. 

********************************************************************** 


thf following outputs arf characteristic of 
A SINGLE transfer function frequency response. 


-L-omtpUT matrix -ar- 

this is the pfouceo a* matrix 

FOR A PARTICULAR USER SPECIFIED 
TRANSFER FUNCTION TYPE . 

PFF. input data list and 

SECTION ITT. 0-2. IN VOL I. 

THF ROOTS OF AP ARE THE 
TRANSFER FUNCTION POLFS. 


OUTPUT MATRIX BCOl 

THIS IS THE VECTOR (COL) WT TH WHTCW 
AP IS AUGMENTFO TO detfRminf THF 
TRANSFER FUNCTION ZEROS. 


* 
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deal part imaginary part rfal PART IMAGINARY part 


ThFSF ARF THF COMPLEX POOT a rpays 
ae: FXTRACTFn FROM MATRICES 
ar and AR TRANSPOSE RESPECT TVFI Y . 

Thf USER SELECTS VIA INPUT WHTCH 
Srj OE ROOTS TO USE FOR THF POIFS. 


NUM 

REAL PART IMAfiTNARY PART 

* 

numerator roots 

(TRANSFER FUNCTION 7ER0S j 


OEN 

RFAL PART IMAGINARY PART 


nENOMTMATOR ROOT* 
(TRANSFEP FUNCTION PO|F«=) 
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c 

c 

r 

c 

c 

c 

c- 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

f 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


OUTPUT MATQIX rpfo 

transfer function poot appay containing BOOT, 

COUNTS* TIMF CONSTANTS* DAMPING* AMD FPFOUFNCY 
FOR ZFROS ANO POl ES. 


FI.F 

1 

E 

NO. 

OF 

NUMFPATOR REAL ROOTS, 

MNP 

fle 

? 

E 

MO. 

OF 

MUMFRATOR COMPLEX PATPS* 

NNC 

FLF 

3 

S 

MO. 

OF 

mijmfpator frfe s*s* 

NN7 

ELE 

4 

= 

MO. 

OF 

nFNOMTMATOR 

RFAL ROOTS. 

NOR 

FLF. 

5 

S 

MO. 

OF 

hfnomtmator 

complfx datrs 

*MDC 

flf 

<S 

= 

MO. 

OF 

OFNOMTMATOR 

FPFF S»S» 

NOZ 

FLE 

7 

x 

POOF GAIN, k-P. 




THF boots FOLLOW in THF OPOFQ 

ELFMFNT LOCATION OFSCRIPTION 

P THRU NNR*7 NUMFRaTOP TI^F CONSTANTS 

NNP*ft THRU ?*NNC ♦ T NUMFPATOR DAMPING AND FPEOUFNC TFS — 

ZFT Alt OMFGAl* 7 FTaP, OMFOA? 

FPFF S»S ARE NOT INCLUDED 

REMAINING FI.FMFNTS APF DFMONINATOR ROOTS TN SAMF 
ORDER AS NUMFPATOR ROOTS, TF,* 

TaU<1) *..*TaU(NDB) fZFTAU) * OMEGA (1) t • • • ZFT A fNDO * OMEGA (NDC) 


« 
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c 

r 

c 

c 

c 

r 

r 

c 

c- 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r 

r 

c 

c 

r 

c 

c 

r 

c- 

c 

c 

c 

c 

r 

c 

c 

r 

c 

c 

c 

r 


ThP transffq FUNCTION 

FRFQ/PAO/SFC FRF0/WERT7 


frfoufncv PFSPOnSF Fallows 
RFAt. I HAG AMP nFCTRH S 


p AD 


ThF DAMPED PFCONANCFF (RftSFD ON a 
POLE OR 7FR0 ) APF IDENTIFIED 

PY ***** ROTH THF LFFT AND RIGHT MaOGTNS. 


****************************************** 


THF FOLLOWING OUTPUTS ARE TYPICAL OF A POOT 
LOCUS INVESTIGATION. 


output matrices POFN 

PNU M 

ThF HFNOMINATOP and NUMERATOR 
transfer FUNCTION POLYNOMTNAI 
COFFFICTFMTS ASCFNnlNG POWERS OF s. 

P(S> 3 NUMFP&TOR 

ANn 

0(G) = 0FMONTNAT0R 

PPFPPOCE SSEO 

POLYNOMIAL cofffictfnts as used py RLOCUS. 


DFG 


C 

= USER TOFNTTfTED STARTING POTNT 
FOR The LOOT. 

= I IMITS On real AND imaginary 
COMPONENTS FOP THF LOCI. 


roots frror 


ROOT LOCUS OUTPUT 

c 

c 

c 

c**»* *************** ******************** ** * **************************** 


c 

C 

r 

c 

c 

c 

c 

c 

c 

c 

c 


STARTING POINT 
SCAN LIMITS 


gain 


ORIGINAL PAGE IS 
OF POOR QUALITYi 


* 
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VI. 


AUXILIARY PROGRAMS 


This section describes two auxiliary digital codes that have 
been developed to aid the DISCOS program system user. The first 
code is a FORTRAN program which accepts the DISCOS code as input 
and, based upon some additional user-supplied input, automatical- 
ly redimensions the DISCOS source program to minimize core stor- 
age requirements. The second code is a DISCOS /NASTRAN interface 
which processes user-supplied NASTRAN generated data into the 
required DISCOS input formats. 


A. REDIM - THE REDIMENSION PROGRAM 


This code was developed to aid the user in the efficient use of 
available digital computer core storage locations. Examination 
of existing digital computer codes for generalized analyses of 
(possibly) large systems indicates that very frequently the 
nature of the code dictates that a great deal of core storage 
locations are required (due to the sizes of program DIMENSION 
and COMMON blocks) . This often leads to inefficient use of 
core storage as the user must have available sufficient core 
storage locations so as to satisfy the program size. As a 
large percentage of program executions probably don't require 
the maximum dimension sizes of program storage blocks, it is 
obvious that a automatic procedure to alter the program code to 
meet a user's specific requirements would be desirable. Program 
REDIM was developed to satisfy these requirements. 

REDIM is a self-contained code that contains an extensive list 
of format statements. The code reads the DISCOS source code 
from tape as coded data and reproduces it on the tape unless 
it finds an identifying format number in columns 73-75. In 
this case, it rewrites the source code according to the format 
corresponding to the identification. REDIM, therefore, provides 
an efficient and foolproof method of recasting the source input 
code to meet the user's requirements. 

Following is a more detailed explanation of the manner in which 
the REDIM program can be used . 
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oooooonoooooon 


***************$ *******^* ***************** **** 5 ** *************** ******* 
**************** ******************************** *********************** 


PROGRAM REDIM INPUT DATA STREAM 


★****♦**• 


THERE is NO CARC INPUT data 0 toil RED fop THIS 
PROGRAM. IF A USER WISHES TO RED I MENS I ON THE 
DISCOS PROGRAM VIA THE REDIMENSION PROGRAM 
HE MUST MAKE THE NECESSARY ADJUSTMENTS TO 
THE VARIABLES INDICATED BELOW <AnD DEFINED 
A PROGRAM DATA STATEMENT), RECOMPILE A\C 
EXECUTE THE RED IM PROGRAM 


BY 


THE FOLLOWING 
ADJUSTED 


PROGRAM DATA VARIALLES CAN 


CURRENT VALVE 


£**$* 9 $$* *# ****#£**** * *** ** ***#*********##* + ***jMc**S**$******#a*##5*:ir*:$t** 
£* *********************** ***** + ******* *********** »#*^**** #**--*#*#»*:?<* 

C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


NBMAX 


MAXIMUM 

NUMBER 

OF 

BODIES 

6 

NHMAX 

= 

MAXIMUM 

NUMBER 

OF 

HINGES 

6 

NSPMAX 

= 

MAX I MUM 

NOKEER 

OF 

SENSOR POINTS 

IS 

NMWMAX 

= 

MAXIMUM 

NUMBER 

OF 

MOMENTUM, WHEELS 

n 

NMWEOD 

— 

MAXIMUM 

NUMBER 

OF 

MOMENTUM 




WHEELS 

PER BODY 


4 

NMDBOD 

ZL 

MAX IMUM 

NUMBER 

OF 

ELASTIC 




MODES 

PER PCD 

Y 


I 2 

KY 

= 

MAXIMUM 

NUMbER 

OF 

STATE VARIABLES 

25 C 

JMAXC 

= 

MAXIMUM 

NUMBER 

OF 

JOINTS LN A 




EDDY - 

CONSIST 

ENT 

MASS 

7 

JMAXL 

— 

MAXIMUM 

NUMBER 

OF 

JOINTS ON A 




BODY - 

LUMPED 

MASS 

105 


» 
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CURRENT valve 


MAXZP 


ROW DIMENSION CF TIME HISTORY 

PLOT DATA ARRAY 1000 


MXJVPL = SIZE OF TIME HISTORY PLOT 

VARIABLE SELECTION VECTOR — MUST 
NOT BE CHANGED — U 


MAXDUM = SIZE OF DUMMY VECTOR IN TIME 
HISTORY PLOT SECTION -- MUST 
NOT b£ LESS ThAN NCPLt T — 1500 


MAXCNT = SIZE CF CONTROL DATA VECTOR 


100 


LI 


L *> 


= ROW DIMENSION OF WORK SPACr 
THAT WILL ACCOMMODATE THE 
LINEARIZED COEFFICIENT MATRIX 
(A) IN DYNSau SEGMENT IOC 

- SIZE CF FFLCJENCY RtSRCNSl' 

DATA VECTORS AND CORK t SPEND S 
TO MAXIMUM NC. OF POINTS TO 
EE PLOTTED 50G 


****** ******************************* ****** ***** **** **** » ** .1 *********** 
0 .*****.* ****************************************** ************ *********** 
c 
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NASFOR - THE NASTRAN INTERFACE PROGRAM 


The multi-purpose programming system described herein can be 
made more versatile if a reliable and efficient means to pro- 
cess input data arising from other sources can be provided. 

One other source of input data is NASTRAN, a digital code that 
is gaining wide acceptance in the aerospace and other industries. 
Program NASFOR, described in this section, was developed to pro- 
vide an interface between NASTRAN and DISCOS. The program as- 
sumes NASTRAN generated structural data is available in a 
prescribed format and transforms these data to a format accep- 
table to the DISCOS system. Either tape or punch card data 
may be processed. 

NASFOR is a self-contained code; it processes data from one 
source (NASTRAN) and generates data for application in DISCOS. 
Originally, it was felt that this interface should be an inte- 
gral part of the DISCOS code. However, during development, it 
was realized that this would impose a large overhead on the 
dynamic response program and it was, therefore, decided that 
NASFOR should be a stand-alone program. 

NASFOR has the capability to process either tape or punch card 
input and create either tape or punch card output. The output 
formats are consistent with the input requirements of the DISCOS 
program system. The input formats assume that NASTRAN generated 
data is double precision and in 0UTPUT2 format if on tape or is 
single precision and of the format (24X.3F8.0) if on punched 
cards. Reference to the example following indicates the format 
requirements for all other date required to exercise the pro- 
gram. 

The code was designed to process NASTRAN generated structural 
data for a series of bodies and to create DISCOS input compati- 
ble with subroutine MSMODL, the lumped mass input routine. It 
is assumed that the available data is in a specific format as 
follows: 

For a body whose inertial and geometric characteristics are 
defined at a set of NJ discrete joints, 



1. Inertial Properties 








m 


s 

z 

-s 

y 



m -S z 


s 

X 



m S y 

- S x 


m i = 


J 

XX 

-J 

xy 

-J 

XZ 



(Syra) 

J 

yy 

-J 

yz 


- 



J 

zz 

where 

m 

= mass 




S 

= static mass moments 


x.y.z 

j * inertias 

XX j i • * j zz 

and the total assembled mass matrix is of the quasi-diagonal 
form 



M = 


m. 


“NJ 


6NJx6NJ ' 


2. 

Joint Coordinate 

Locations 


i) 

card input data 



G = 

[x y 

Z NJx3 



with x = 

J 

X 1 x 2 ' 

• • • *NjJ 

T 

lxNJ 


y = 

K y 2 • 


T 

lxNJ 
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“2 ' ' ' 
or ii) tape input data 

G * G. - G 
1 o 

„ith Gj - In x y x NJx4 

X 

where n ■ |_ni n 2 . . . n N jJ ]_ xN j are joint numbers and will be 
ignored by tne tape reading section and x, y and z are as above 

and G * x y z 

L J NJx3 

where x Q , y Q and z Q are user-supplied card inputs. and may be 
null. 

3. Modal Properties 

* “ * R *eI 

L J 6NJxNM 

where $> R = NR rigid body modes 

4> ■ NE elastic modes and NM = NR + NE 

£ 

and where 

*j ■ [*i *2 ; ■ • *i • • • *>«J i x61u 

and * t - |^h x h y h z „ y „ 2 J 

with h “ modal displacement amplitude 

x,y,z 

a - modal rotation amplitude. 

x,y,z 
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4. Generalized Stiffness and Damping 


K = K and C 

gen J NExNE 



NExNE 


The data input as previously noted, are manipulated within the 
program and the results are written on tape and/or provided as 
punch card output as follows! 


1. inertial properties 





NJx3 


J = 


xx ^ 


L J xxNJ 


J yyi 

J yyNj 


J zz\ ^xy^ J xz^ 


J zznj J xyNJ J xzNJ 


J yzi 


J yzNJj 


NJx6 


2 . joint coordinate locations 


G = 


X 1 y l Z 1 


X NJ y NJ Z NJ. 


NJx3 



3. modal properties 


x,y,z 


h 1,1 h 1,2 

x,y,z x,y,z 


, h ^’ 1 

L x,y,z 


. h 


1,NE -| 
x,y,z 


NJ ,NE 
x,y,z J 


x,y,z 


iV.i 


u x,y ,z 


1,2 


x,y,z x,y,z 


NJ , 1 

T 7 


1 ,NE 


• ^x.y.z 


NJ ,NE 
x.y.z 


4, generalized stiffness and damping 


K 


gen 


K 11 K 12 K 1,NE 


^,1 • ' • 1 ^NE,NE 


J NExNE 


gen 


C 11 C 12 


"'NE , 1 


'1,NE 


J NE,NE 


J NExNE 


NJxNE 


NJxNE 



NASFOR consists of a main program, 8 program subroutines to pro 
cess the data and 10 auxiliary input/output routines. The func 
tion of the program routines is indicated in Table VI.B-1. and 
logic flow diagram appears as Figure VI.B-1. 


Table VI.B-1. Description of NASFOR Subroutines 


Subroutine 

Func tion 

MAIN 

TFETCH 

CFETCH 

GMASS 

GMODE 

GGEGM 

GSTIF 

GDAKP 

PRINTT 

Program control 

Fetch a matrix from NASTRAN tape 

Fetch a matrix from cards 

Generate mass data 

Generate modal displacement data 

Generate geometric data 

Generate stiffness data 

Generate damping data 

Print entire input tape (on option) 


The code has been designed to minimize input data requirements 
yet provide a high degree of flexibility. Following is a de- 
tailed explanation of the input data stream requirements. 





999V CALL START 


CALL PAGEHD 


WRITE {NOT, 1) 


CALL C GHENT 


READ (NIT, 1000) NTAPE1, NTAPE2, NTAPE3 
READ (NIT, 1001) TAPEID 

READ (NIT, 1000) IFPRNT1, IFPRNT2, IFPRNT3 , IFPNCH 


WRITE (NOT, 1002) NTAPE1, NTAPE2, NTAPE3, TAPEID 
WRITE (NOT, 1003) IPPRNT1, IFPRNT2, IPPRNT3, IFPNCH 



WRITE (NOT, 991) 


NT APE 2 - LAIS (NTAPE2) 


REWIND NTAPE1 


READ (NIT, 1000) NB 


DO 5000 NBODY - 1, NB 


CALL PACEHD 


READ (NIT, 1000) NJ , NE, NR 


WRITE (NOT, 

1004) 

NBODY, NJ, NE, NR 

READ (NIT, 

1003) 

( (IDMAT (I, J), J-l, 5), 1-1, 2) 



n 

■ i 

o o 


IDMAT (2, Jr 
V - TAPE ^ 






4 CONTINUE 














































MATID 














c 

c 

c 

c 

c 


o*** ****************************************************************** 

c**** ****************************************************************** 

c 

c 

C PROGRAM NASFOR INPUT DATA STREAM 

C 

C 

C ******* 

c 


************************************************************ ******* 
. c**************** ************************************************ ******* 

c 

c 

c 

C PROGRAM READS NASTRAN OUTPUT TAPE AND/OR CAROS, 

C MANIPULATES THE DATA AND WRITES (PUNCHES) TAPE 

C (CARDS) TO BE USED AS INPUT TO PROGRAM DISCOS 

C 
C 

C PURPOSE — PROCESS NASTRAN STRUCTURAL DATA 

C FOR A SERIES OF FLEXIBLE BODIES 

C AND GENERATE DATA CONSISTENT 

C WITH DISCOS SUBROUTINE MSMODL 

C INPUT REQUIREMENTS 

C 
C 

C ASSUMES — FOLLOWING NASTRAN DATA AVAIL ABE 

C FCR ALL FLEXIBLE BODIES 

C 


c 

c 

MATRIX 

SUE 

c 

c 

MASS - LUMPED 

6NJ X 6NJ 

c 

MODES 

6NJ X(NE+NR) 

c 

GEOMETRY 

NJ X A 

c 

STIFFNESS - MODAL 

NE X NE 

c 

DAMPING - MODAL 

NE X NE 


C 

C 
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WHERE NJ- = NO. OF JOINTS ON BODY 

NE = NO. OF ELASTIC MODES RETAINED FOR BODY 
NR =■ NO. OF RIGID BODY MODES 


ASSUMES — 

1. NASTRAN MATRICES WRITTEN DOUBLE PRECISION 
{ RE AL*8 OR I TYPE =2) IF ON TAPE 

2. NASTRAN MATRICES WRITTEN SINGLE PRECISION 
(FORMAT = 24X,3F8.0) IF ON CARDS 


3. MAX(NJ) = 1000 

MAX (NE ) =■ 10 — REDIMENSION IF SO REQUIRED 

MAX (NR) = 6 


NOTE - 

C INDICATES WHERE DATA INPUT 

IS REQUIRED 

NIT = INPUT TAPE NUMBER 


9999 CALL START (SEE FOLLOWING) 


CALL COMENT (SEE FOLLOWING) 


READ (NIT, FORMAT « 315) NTAPE1, NTAPE2* NTAPE3 


NTAPE1 * LOGICAL UNIT NASTRAN INPUT TAPE (EG, 20) 
NTAPE2 * LOGICAL UNIT OUTPUT TAPE (EG, 30) 
NTAPE3 ' SCRATCH (EG, 40) 


- READ(NIT, FORMAT = A6) TAPE ID 

TAPE ID = 6 CHARACTER TAPE IDENTIFICATION (EG, L12345) 
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READ(NIT, FORMAT * 415) IFPRT1 * IFPRT2* 


IF PRT3 t IFPNCH 


IFPRT1 NE 0 
IFPRT2 NE 0 
IFPRT3 NE 0 
IFPNCH NE 0 


PRINT NASTRAN DATA 
PRINT OUTPUT DATA 

PRINT NASTRAN HEX TAPE IF ERROR FOUND 
PUNCH OUTPUT DATA 


READ(NIT, FORMAT = 15) Nb 

NP = NO. OF BODIES FOR WHICH DATA AVAILABLE 
ON THE NASTRAN TAPE (EG, 2) 


I F (NE .EQ. 0) STOP WRITE ERROR MESSAGE 

AND TERMINATE 


THE FOLLOWING LOOP IS EXECUTED FOR EACH BODY 


DO 5000 NBCDY = 1,NB 


READ (NIT, FORMAT = 315) N J , NE , NR 

NJ =• NO. OF JOINTS ON BODY (EG, 11) 

NE = NO. OF ELASTIC MOOE RETAINED FOR BODY (EG, 5) 
NR - NO. OF RIGID BODY MODES (EG, 0) 


— — READ(N1T, FORMAT = 5(A4,6X) ) ( ( IDMAT ( I , J ) , J = 1 , 5 ) , I = 1 ,2 ) 

IDMAT IS AN INPUT ARRAY OF SIZE 2X5 THAT 
DEFINES THE ORDER OF INPUT DATA TO BE READ 
FROM TAPE OR CARDS BY MATRIX FUNCTION AND 
ALSO SPECIFIES WHETHER DATA IS TO BE READ 
FROM TAPE OR CARDS 

THE DATA WILL BE PROCESSED IN THE ORDER J=l,2,..*5 
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AN EXAMPLE IDMAT ARRAY IS 


*** 




*** 

IDMAT = * 




* 

*GEOM 

MASS 

MODE 

STIF 

DAMP* 

* 




* 

♦CARD 

TAPE 

TAPE 

TAPE 

TAPE* 

♦ 




* 

*** 




*** 


AND INDICATES DATA TO BE PROCESSED AS 

1. JOINT GEOMETRIC DEFINITION FROM CARDS 

2. JOINT LUMPED MASS DATA FROM TAPE 

3. ELASTIC MODAL DATA FROM TAPE 
A. MODAL STIFFNESS FROM TAPE 

5. MODAL DAMPING FROM TAPE 


THE FOLLOWING LOOP IS EXECUTED 5 TIMES PER BODY 


DO 4999 MAT = 1*5 
MATFMT = IDMAT ( 2 *MAT) 

IF (MATFMT .EQ. TAPE) CALL TFETCHC ) 

IF (MATFMT .EQ. CARD) CALL CFETCH 

TFETCH IS TAPE PROCESSING CONTROL SUBROUTINE 
CFETCH IS CARD PROCESSING CONTROL SUBROUTINE 

SEE FOLLOWING FOR DESCRIPTION 

4999 CONTINUE 

5000 CONTINUE 
GO TO 9999 
END 
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*********************************************************************** 

* 

* SUE ROUT I NE START - INPUT IDENTIFICATIONS 

* 

****** ************ ******************************************** ********* 

READINIT , FORMAT = A6,4X,3Afc) 1RUNNO (UNAME (1) t 1-1 .3 ) 

I RUN NO = RUN IDENTIFICATION (6 CHARACTERS) 

UNAMF = USERS NAME (16 CHARACTERS) 

IRUNNO EQ 4HST0P - TERMINATE THE RUN 
IRUNNO NE ‘rHSTOP - CONTINUE THE RUN 


READ(NIT, FORMAT » 12Ao) (TITLE1 ( 1 ) » 1-1 » 12 ) 

READ(N1T, FORMAT = 12A6) (TITLE2 ( X) , 1 = 1 , 12 ) 

TI7LE1 = 72 CHARACTER TITLE 
TITLE2 = 72 CHARACTER TITLE 

RETURN 

END 


*********************************************************************** 

* 

* SUBROUTINE COMENT - INPUT USER SUPPLIED COMMENTS 

* 

*********************************************************************** 

READ(NIT, FORMAT = 13At>,lX,Al) ( IREMRK 1 1 ) , 1*1 , 13> , IPGHD 

IREMRK * 76 CHARACTER COMMENT 
IPGHD * NEW PAGE FLAG 

IPGHD * 1 HP — NEW PAGE BEFORE PRINTING 
THERE IS NO LIMIT TO THE NUMBER OF COMMENT 
CARDS BUT THE LAST ONE MUST CONTAIN ZERO 
IN COLUMNS 1 THRU lo 

RETURN 

END 
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C 

£****** ************************************* ********************* ******* 
c* 

C* SUBROUTINE TFECTH 

C* ********************************************* 

THIS SUBROUTINE READS NO INPUT DATA 
DIRECTLY BUT DOES CALL AUXILIARY 
SUBROUTINES DEPENDING UPON THE VALUE 
OF THE INPUT ARGUMENT -MAT1D- 

THESE AUXILIARY SUBROUTINES READ 
DATA AS DESCRIBED LATER 

IF (MATID *EQ. 1) CALL GMASS ( . . ) LUMPED MASS DATA 

IF (MATID .EQ. 2) CALL GMODE ( • • ) MODAL DATA 

IF (MATID .EQ. 3) CALL GGEClMU.) GEOMETRIC DATA 

— ———————— IF (MATID • EQ • A) CALL GSTIF(.«) STIFFNESS DATA 

IF (MATID .EC. 3) CALL GDAMP ( . • > DAMPING DATA 

RETURN 
END 


SUBROUTINE CFETCH 

************************************************ *********************** 

THIS SUBROUTINE READS CARD INPUT DATA AS FOLLOWS — 

READ (NIT* FORMAT = A6,4X,215) ANAME, NR, NC 

ANAME * 6 CHARACTER MATRIX IDENTIFICATION 
NR * NO. OF ROWS IN MATRIX 

NC * NO. OF COLS IN MATRIX 

DO 100 I = 1 »NR 

REA0(NIT, FORMAT * 24X,3F8.0) (WR( J ) , J=1 ,NC ) 

( WR ( J )«J S 1*2,... * NC) * INPUT ROW I OF MATRIX 

100 CONTINUE 

RETURN 
END 


VI-18 



c 
c 
c 
c 
c 

£****♦***********: 

o 

c* SUBROUTINE GMASS 

C* 

C* ****************************************** ******** ******** ******* 

c 

c SUBROUTINE READS OUTPUT 

EACH BODYS MASS, STATIC 
INERTIAL DATA 


TITLE CARD FOR 
MASS MOMENT AND 


DATA PEAD ON FIRST CALL 
RE AD (NIT , FUKM AT = 3(A6,AX)) 


FOR EACH BODY 
AN1 * AN2, AN3 


AN! 
AN 2 
AN3 

RETURN 

END 


CHARACTER OUTPUT TITLE 

character output title 

CHARACTER OUTPUT TITLE 


FOR MASSES 

FOR STATIC MOMENTS 

FOR INERTIAS 


C 
C 
C 
C 
C 

C 

C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

c****** ******** ********************************************** *********** 

c* 

o SUE ROUTINE GMCDE 

C* 

c ************************************** * ***************** *************** 

c 

SUBROUTINE READS OUTPUT TITLE CARD FOR 
EACH BODYS MODAL DATA 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


DATA READ ON FIRST CALL 
READ(N1T, FORMAT * e> ( A r* » AX ) ) 


FOR EACH BODY 

AN1, AN 2, AN3 , ANA, AN5 , 


AN6 


AN1 
AN* 
AN 2 
ANA 

AM 
AN 6 


CHARACTER 
CHARACTER 
Ch AR ACT r R 
CHARACTER 
CHARACTER 
CHARACTER 


CO T FUT TITLE 
OUT 'UT TITLE 


CUTFc T 
OUTPUT 
OUTPUT 
OUTPUT 


TITLE 

TITLE 

TITLE 

T a TLE 


RETURN 

END 


FOR hX MODAL AMPS 
FOR HY MODAL AMPS 
FOR hZ MODAL AMPS 
FOR SIGX MODAL SLOPES 
FOR SIGY MODAL 
FOR SIGZ MOOAL 
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SLOPES 

SLOPES 



C*********************************************************************** 


C* 

c* SUBROUTINE GGEOM 

C* 

C**************** ******************************************** *********** 


c 

c SUBROUTINE READS INITIAL GEOMETRIC DATA TO 

C BE MANIPULATED WITH INPUT TAPE GEOMETRIC 

C DATA AND OUTPUT TITLE CARD FOR EACH 

C BODIES GEOMETRIC DATA 

C DATA READ ON FIRST CALL FOR EACH BODY 

C INPUT FORMAT FOR SUBROUTINE READ EXPLAINED 

C ELSEWHERE IN THIS DOCUMENT 

C 

c INITIAL GEOMETRIC DATA (MAY BE NULL) IS 

C PROVIDE TO ACCOUNT FOR POSSIBLE 

C GEOMETRIC OFFSETS 

C 

c CALL READ ( WS ,NR ,NC , KWS , LWS ) 


C 

C 

C 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c- 

c 

c 

c 

c 

c 

c 


NR = INPUT NO. OF ROWS (MUST = NJ) 

NC = INPUT NO. OF COLS (MUST = 3 ) 

WS = INPUT MATRIX OF INITIAL 

GEOMETRIC COORDINATES — SIZE = NJ X 3 


THE FINAL GEOMETRY FOR NJ JOINTS 
WILL BE COMPUTED AS 


*** 

*** 

*** 

*** 

*** 

*** 

* X ( 1 1 

YU) 2(11 ^ 

* 

* 

* 

* 

* . 

• . * 

* X Y 

Z * 

♦ X 

Y Z * 

* . 

• . * 

=. * 

♦ - 

♦ 

* 

* . 

• * * 

* (TAPE 

) * 

* (CARD) ♦ 

* . 

# , * 

* 

* 

* 

* 

* . 

• • * 

* 

* 

* 

* 

* X ( NJ ) 

Y( NJ ) Z(NJ>* 

♦ 

* 

♦ 

♦ 

*** 

*** 

*** 

*♦* 

*** 

**** 

SO 

THAT, FOR THE ITH 
XU) * X(I) FROM 

JOINT 

TAPE INPUT 

- X(I) 

READ 

ABOVE 

(EAD(NIT 
AN1 = 

Y ( I ) = Y ( I ) FROM TAPE INPUT 
Z(I) = Z(I> FROM TAPE INPUT 

t FORMAT * A6) AN1 

6 CHARACTER OUTPUT TITLE FOR 

- Y(I) READ ABOVE 

- Z(I) READ ABOVE 

GEOMETRIC DATA 


RETURN 

END 
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♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦I************************************************ 

* 

* SUBROUTINE GST1F 

* 

************* ***************** **************************** ****** ******* 


SUBROUTINE READS OUTPUT TITLE CARD FOR 
EACH BODYS MODAL STIFFNESS DATA 

DATA READ CN FIRST CALL FOR EACH BODY 


READ (NIT, FORMAT = A6 ) AM 

AN 1 * 6 CHARACTER OUTPUT TITLE FOR 
MODAL STIFFNESS DATA 


RETURN 

END 


***** ** ******** ************************************ ******** *********** 

* 

* SUBROUTINE GDAMP 

* 

*********************************************************************** 


SUBROUTINE READS OUTPUT TITLE CARD FOR 
EACH BODYS MODAL DAMPING DATA 

DATA READ CN FIRST CALL FOR EACH BODY 


READ(N1T, FORMAT * A6 ) AN1 

AN1 = 6 CHARACTER OUTPUT TITLE FOR 
MODAL DAMPING DATA 


RETURN 

END 
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APPENDIX A — INPUT/OUTPUT SUBROUTINE EXPLANATIONS 


This appendix presents excerpts from Synthesis of Dynamo Sys- 
tems Using FORMA — Fortran Matrix Analysis, MCR-71-75, Martin 
Marietta Corporation, Denver, Colorado, May 1971, that explain 
the subroutines from the FORMA library used in the digital com- 
puter program. 


A-l 


COMENT 


Subroutine COMENT reads input comment cards and reproduces 
each card in the printed output of the computer run. Each com- 
ment card may have any keypunch symbol in card columns 1 thru 
78. A use of COMENT is to print an explanation of coordinates 
used in a computer run. Thus, this information is always re- 
tained with a run to correlate matrix location numbers with 
physical coordinates. 


R, L. Wohlen 
May 1971 

Last Revision: R. L. Wohlen 

Sept. 1971 
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INTAPE 


Subroutine INTAPE initializes a tape (a disk is preferred. See 
writeup of Subroutine WTAPE. ) for the FORMA tape system by writing 
EOT (end of tape) at the beginning of the tape (disk). All FORMA 
tape subroutines recognize this EOT as being the end of written 
data. Each "new" tape (disk) must be initialized with this Sub- 
routine INTAPE to make the tape (disk) compatible with the other 
FORMA tape subroutines (LTAPE, RTAPE, WTAPE, and UPDATE). 

A "new" tape (disk) is defined as a tape (disk) for which it 
is desired to s'.art writing matrix data at the front of the tape 
(disk). Thus, a "new" tape (disk) could be one with obsolete 
FORMA matrix data on it as well as one that has never been written 
on by the FORMA system. 

As an example, pertinent statements from a program containing 
INTAPE could be: 

DATA NIT, NOT/ 5, 6/ 

1001 FORMAT (12A6) 

NRTAPE = 10 

READ (NIT, 1001) IFINIT, TAPEID 

IF (IFINIT .EQ. 6HINITIL) CALL INTAPE (NRTAPE, TAPEID) 


The input data (starting in card column 1) to this example 
program would be: 

either INITILTXXXX, if the tape is to be initialized. 

(TXXXX represents the particular tape number 
used, e.g., T1234) ; 

or NOINIT, if the tape is not to be initialized. 

The tape identification is not needed. 


R. L. Wohlen 
May 1971 



LTAPE 


Subroutine LTAPE lists the matrix headings (see Subroutine 
WTAPE writeup) written on a FORMA tape (or disk). These matrix 
headings were written by Subroutine WTAPE and consist of: 


NO. 

RUN NO. 

NAME 

NROWS 

NCOLS 

DATE 

NNZ 


= Matrix number on tape; 

= Run number of problem when matrix was written on 
tape ; 

= Matrix name; 

* Number of rows of matrix; 

* Number of columns of matrix; 

■ Date when matrix was written on tape; 

■ Number of nonzeros (just used in sparse FORMA 
where only nonzeros are used) ; 


PARTITION = Partition number of sparse matrix. 
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READ— 1/6 


Subroutine READ reads a matrix of real numbers (a FORTRAN 
term for numbers with a decimal point) from either cards or tape 
into the computer. The matrix is then printed so that these in- 
put data are recorded with the answers of a run. A print sup- 
pression option is available for a matrix read from tape. On 
option, the matrix read from either cards or tape may be written 
on a tape (by Subroutine WTAPE). 

The first data card read by Subroutine READ contains the in- 
formation to indicate whether cards or tape will be used. The 
information entered on this card (and subsequent cards for card 
input) is given below. 


Card Data Input Form 

Required entries are denoted by an * symbol below. Any other 
entry is optional 


First Card 


Card 

Columns 

1-6 

7-10 

11-15 

16-69 

72 


Format 
Type (1) 


En t vy 


(0 

§ 

•H 

4J 

a 

o 

0) 

8* 

V 

0) 

u 

■H 

U 

3 


72 


73-78 


7 3-78 


73-76 


73-78 

79-80 


A ^Matrix Name. Will appear 

in printout. 

I ^Matrix Row Size. 

I ^Matrix Column Size. 

A Any remarks to further iden- 

tify the input matrix. 

$. Only if the Write-Tape 
is to be initialized by Sub- 
routine INTAPE. The Write- 
Tape identification will be 
* from card columns 73-78. 
or Anything other chan $ is the 

Write-Tape is not to be ini- 
tialized. 

A The Write-Tape identification. 

(e.g. , T1234) . Use with $ in 
card column 72 . 

or REWIND. The Write-Tape will 

be rewound before being used, 
or LIST. The Write-Tape will be 

listed by Subroutine LTAPE 
after the matrix has been 
written on the Write-Tape, 
or Anything else will be ignored, 

1 The Write-Tape Number, (e.g., 

or Blank if the matrix is not to 
be written on tape. 
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Card 

Format 


Columns 

Type (1) 

Entry 

Middle Cards 1-5 

I 

*Row Number of matrix elements 
on card. 

6-10 

I 

^Column Number of matrix ele- 
ment in first data field. 

11-27 

E 

*First data field with matrix 
elements. (2) 

28-44 

E 

^Second data field with matrix 
elements. (2) 

45-61 

E 

*Third data field with matrix 
elements. (2) 

62-78 

E 

^Fourth data field with matrix 
elements. (2) 

Last Card 1-10 

I 

*Ten zeroes. 


Note (1) Format Type A allows any keypunch symbol. 

Format Type I allows only integer numbers right justified 
in the field. Format Type E allows only real numbers 
(a FORTRAN term for numbers with a decimal point) any- 
where in the field. 

Note (2) Only nonzero elements need be entered. 


As an example of card input to Subroutine READ consider the 
following matrix: 


[ a1 *c] 3x6 


1. 0. 3. 0. 6. 

0. 2. 4. 0. 0. 

0. 7. 0. 0. 0. 



0 . 

0 . 

mm 


This matrix is also to be written on tape number 21 that is to 
be initialized and identified as T4334. Figure 1 demonstrates 
how this information could be written on a coding form to facili- 
tate keypunching to cards. 
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Figure 1 Example of Card Input for Subroutine READ 
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Tape Data Input Form 

Required entries are 

denoted with an * symbol below. Any 

other entry is optional. 

Only one 

card is used for each matrix 

read. 

Card 

Format 


Columns 

Type (1) 

Entry 

One Card 1-6 

A 

*Name of matrix to be read 
from the Re ad -Tape. 

10 


Zero. The Read-Tape will 



move forward from its present 
position and search to the 
. end of the tape. If the 
matrix is not found upon the 
first end-of-tape encounter, 
the tape will automatically 
rewind and make one more 
pass. If it is not found on 
the second end-of-tape en- 
counter, an error message 
will be printed and the pro- 
gram will stop. 

7-10 

I or 

Minus the location number of 
matrix on the Read-Tape. Tap< 
will be positioned at the be- 
ginning of the location sped 
fied and then continue as 
described above for a zero 
in column 10. 

11-15 

I 

*The Read-Tape Number, (e.g., 
If positive, the matrix read 
will be printed in the output 
If negative, the matrix read 
will not be printed in the 
output. 

16-21 

A 

*Run number of matrix to be 
read from the Read-Tape. 

22-27 


REWIND. The Read-Tape will 
be rewound before being used. 

22-25 

or 

LIST. The Read-Tape will be 
listed by SuDroutine LTAPE. 

22-27 

or 

Anything else will be consid- 
ered as part of the remarks 



described below. 
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Entry 

Any remarks to further iden- 
tify the input matrix. 

$. Only if the Write-Tape 
is to be initialized by Sub- 
routine INTAPE. The Write- 
Tape identification will be 
from card columns 73-78. 

Anything other than $ if the 
Write-Tape is not to be ini- 
tialized . 

The Write-Tape identification, 
(e.g., T1234) . Use with $ in 
card column 72. 

REWIND. The Write-Tape will 
be rewound before being used. 

LIST. The Write-Tape will 
be listed by Subroutine LTAPE 
after the matrix has been 
written on the Write-Tape. 

Anything else will be ignored. 

The Write-Tape Number, (e.g., 21). 
Blank if the matrix is not 
to be written on tape. 

Note (1) Format Type A allows any keypunch symbol. 

Format Type I allows only integer numbers right justi- 
fied in the field. 

As examples of tape input to Subroutine Read consider: 

Example 1. A matrix named AB2 with run number of RUN-46 is to 
be read from tape number 11 into the computer and 
printed. This matrix is also to be written on tape 
number 22 that is to be initialized and identified 
a9 T4321. 

Example 2. A matrix named XYZ4 with run number of TKD is on tape 
number 13 twice. The first time is at location 29 
and the second time is at location 54. It is desired 
to read the second matrix. 

Figure 2 demonstrates how these two examples would be wtitten 
on a coding form to facilitate keypunching to cards. 
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Subroutine READIM reads a matrix of integer numbers from 
either cards or tape into the computer. The matrix is then printed 
so that these input data are recorded with the answers of a run. 

A print suppression option is available for a matrix read from 
tape. On option, the matrix read from either cards or tape may 
be written on a tape (by Subroutine WTAPE) . 

The first data card read by Subroutine READIM contains the 
information to indicate whether cards or tape will be used. The 
information entered on this card (and subsequent cards for card 
input) is given below. 


Card Data Input Form 

Required entries are denoted by an * symbol below. Any other 
entry is optional. 


Card Format 

Columns Type (1) Entry 


First Card 1-6 

7-10 

11-15 

16-69 

r 72 


S 

0 > 

s- 

H 


U 

3 


72 


73-78 


73-78 

73-76 


73-78 

79-80 


A 

I 

I 

A 


or 


A 


or 

or 


or 

I 

or 


♦Matrix Name. Will appear 
in printout. 

♦Matrix Row Size. 

♦Matrix Column Size. 

Any remarks to further iden- 
tify the input matrix. 

$. Only if the Write-Tape 
is to be initialized by Sub- 
routine INTAPE. The Write- 
Tape identification will be 
from card columns 73—78. 

Anything other than $ if the 
Write-Tape is not to be ini- 
tialized. 

The Write-Tape identification, 
(e.g., T1234) . Use with $ 
card column 72. 

REWIND. The Write-Tape will 
be rewound before being used. 

LIST. The Write-Tape will 
be listed by Subroutine LTAPE 
after the matrix has been 
written on the Write-Tape. 

Anything else will be ignored* 

The Write-Tape Number, (e.g., 21). 
Blank if the matrix is not 
to be written on tape. 
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Card 

Columns 

Middle Cards 1-5 

6-10 

11-15 

16-20 

etc 

76-80 


Last Card 


1-10 


Format 

Type (1) E " tr y 

I *Row Number of matrix elements 

on card. 

I *Column N umb er of matrix ele- 

ment ip first data field. 

X *First data field with matrix 

elements. (2) 

I ^Second data field with matrix 

elements. (2) 

I *Fourteenth data field with 

matrix elements. (2) 

I *Ten zeroes. 


Note (1) Format Type A allows any keypunch symbol. 

Format Type I allows only integer numbers right justi 

fied in the field. 

Note (2) Only nonzero elements need be entered. 

As an example of card input to Subroutine READIM consider the 
following matrix: 


[A!*C] 3xg 


1 0 
0 2 
0 7 


3 0 6 5 

4 0 0 0 

0 0 0 0 


This matrix is also to be written on tape number 21 that is to 
be initialized and identified as T4334. Figure 1 demonstrates 
how this information could be written on a coding form to facili 
tate keypunching to cards. 


A-12 










READIM-- 4/6 


Tape Data Input Form 

Required entries are denoted with an * symbol below. Any 
other entry is optional. Only one card is used for each matrix 
read. 


Card Format 

Columns Type (1) 

One Card 1-6 A 

10 


7-10 I or 


11-15 I 


16-21 

A 

22-27 


22-25 

or 

22-27 

or 


Entry 

*Name of matrix to be read 
from the Read-Tape. 

Zero. The Read-Tape will 
move forward from its present 
position and search to the 
end of the tape. If the 
matrix is not found upon the 
first end-of-tape encounter, 
the tape will automatically 
rewind and make one more pass. 

If it is not found on the 
second end-of-tape encounter, 
an error message will be 
printed and the program will 
stop . 

Minus the location number of 
matrix on the Read-Tape. 

Tape will be positioned at 
the beginning of the location 
specified and then continue 
as described above for a zero 
in column 10. 

*The Read-Tape Number, (e.g., 11). 
If positive, the matrix read 
will be printed in the output. 

If negative, the matrix read 
will not be printed in the 
output . 

*Run number of matrix to be 
read from the Read-Tape. 

REWIND. The Read-Tape will 
be rewound before being used. 
LIST. The Read-Tape will be 
listed by Subroutine LTAPE. 
Anything else will be con- 
sidered as part of the remarks 
described below. 
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Card Format 

Columns Type (1) Entry 


28-69 

72 


3 


72 


p. 

o 

V 

§* 

41 

U 

fl 

u 

3 


73-78 

73-78 

73-76 


73-78 

79-80 


A Any remarks to further iden- 

tify the input matrix. 

$. Only if the Write-Tape 
is to be initialized by Sub- 
routine INTAPE. The Write- 
Tape identification will be 
from card columns 73-78. 
or Anything other than $ if the 
Write-Tape is not to be ini- 
tialized. 

A The Write-Tape identification. 

(e.g.. T1234). Use with $ in 
card column 72. 

or REWIND. The Write-Tape will 

be rewound before being used, 
or LIST. The Write-Tape will 

be listed by Subroutine LTAPE 
after the matrix has been 
written on the Write-Tape, 
or Anything else will be ignored. 

I The Write-Tape Number, (e.g., 

or Blank if the matrix is not 

to be written on tape. 


Note (1) Format Type A allows 
Format Type I allows 
fied in the field. 


any keypunch symbol. 

only integer numbers right justi- 


As examples of tape input to Subroutine READIM consider 


Examp le 1 . 


Example 2. 


A matrix named AB2 with run number of RUN-46 is to 
be read from tape number 11 into the computer and 
printed. This matrix is also to be written on tape 
number 22 that is to be initialized and identified 
as T4321. 

A matrix named XYZ4 with run number of TKD is on tape 
number 13 twice. The first time is at location 29 and 
the second time is at location 54. It is desired to 
read the second matrix. 


Figure 2 demonstrates how these two examples would be written 
on a coding form to facilitate keypunching to cards. 


21 ). 
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Figure 2 Examples of Tape Input for Subroutine READIM 








RTAl’K 


Subroutine RTAPE reads a selected matrix frorp tape (disk) into 
the computer core. The matrix to be selected is identified by the 
desired run number and matrix name. This procedure is accomplished 
by searching the matrix headings (see Subroutine WTAPE writeup) 
until a match with the desired run number and matrix name is ob- 
tained and then reading the matrix elements from tape (disk) into 
the computer core. The search starts from the current position 
(does not rewind) of the tape (disk) and proceeds to the EOT (end 
of tape defined in Subroutine WTAPE writeup) . If the desired 
matrix was not found upon reaching the EOT, a rewind is performed 
and one more search to the EOT is made. If the desired matrix is 
again not found, (1) an error message is printed, (2) a listing 
of the matrix headings is printed (see Subroutine LTAPF. writeup), 
and (3) transfer is made to Subroutine ZZBOMB where the program 
is terminated. 
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START 


Subroutine START performs the following operations : 

1) Reads Input Card 1 for the run number (any keypunch 
symbol in card columns 1 thru 6) and the user’s name 
(any keypunch symbol in card columns 11 thru 28) . 

If the run number is equal to STOP (Card columns 1 
thru 4), the run is terminated. 

If the run number is not equal to STOP, the run con- 
tinues In Subroutine START as follows. 

2) Reads Input Card 2 for Title Card i. Any keypunch 
symbols may be used in Card columns 1 thru 72. 

3) Reads Input Card 3 for Title Card 2. Any keypunch 
symbols may be used in Card columns 1 thru 72. 

4) Initializes page number as zero for use In Subroutine 
PAGEHD . 

5) Interrogates computer for the date. 

Run number, date, page number, user’s name, Title Card 1, and 
Title Card 2 are transferred by a COMMON block labeled. L START for 
use in other subroutines PAGEHD, PL0T1, PL0T2 , PL0T3 , and WTAPE . 

Subroutine START is used to start each computer run in the 
FORMA system and will normally be the first subroutine called in 
a computer program. As an example, pertinent statements from a 
program using START could be: 

1 CALL START 


GO TO 1 
END 


R. L. Wohlen 
May 1971 
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WRITE 


Subroutine WRITE writes a matrix of real numbers (a .Fortran 
term for numbers with a decimal point) on paper. A group of up 
to ten consecutive elements from a row of the matrix are printed 
on each line. If all of the elements of a group are zero, printing 
of this line is suppressed. 

Each matrix printed begins on a new page. On each page of 
printout is the page heading given by Subroutine PAGEHD , the name 
o f the matrix, and the^ row size, and column size of the matrix. 

This is followed by the matrix data. On any line of matrix data 
the first integer number is the row number of the matrix elements 
on that line. The second integer number is the column number of 
the matrix element in the first data field. The next group of 
real numbers (up to ten) are the values of the matrix elements. 

This group of matrix elements is given in consecutive column order. 
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Subroutine WRITIM writes a matrix of integer numbers on paper. 
A group of up to twenty consecutive elements from a row of the 
matrix are. printed on each line. If all of the elements of a 
group are zero, printing of this line is suppressed. 

Each matrix printed begins on a new page. On each page of 
printout is the page heading given by Subroutine PAGEHD, the name 
of the matrix, and the row size and column size of the matrix. 
This is followed by the matrix data. On any line of matrix data 
the first integer number is the row number of the matrix elements 
oh that line. The second integer number is the column number of 
the matrix element in the first data field. The next group of 
Integer numbers (up to twenty) are the values of the matrix ele- 
ments* This group of matrix elements is given in consecutive 
column order. 
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I 


Subroutine WfAPE writes matrix data at the end of existing 
written matrix data on a FORMA tape (disk is preferred, see below). 
Each set of matrix data consists of two logical records. The first 
record contains the matrix heading (tape identification, location 
number, run number, matrix name, number of rows of matrix, number 
of columns of matrix, date, and the word "dense ). The second 
record consists of the matrix elements. 

A schematic representation of the tape (disk) is given by the 
following sketch. 


Beginning 

of 

tape (disk) 



where 

- Matrix heading of the i^ written matrix, 

E = Matrix elements of the i written matrix, 
i 

EOT = End of Tape. Data written by Subroutine WTAPE or 

INTAPE that all FORMA tape subroutines recognize as 
being the end of written data. 

Each vertical line is an end of logical record put on by com- 
puter system's routines. The tape is written in binary form as 
opposed to binary coded decimal (BCD) form. 

To find the end of written matrix data, a search is made of 
the matrix headings until the EOT is found. For this reason, a 
"new" tape (disk) must be initialized with Subroutine INTAPE so 
that the tape (disk) contains an EOT. A "new" tape (disk) is de- 
fined to be a tape (disk) for which it is desired to start^writing 
matrix data at the front of the tape. (disk). Thus, a new tape 
(disk) could be one with obsolete FORMA matrix data on it as well 
as one that has never been written on by the FORMA system. When 
the EOT is found, a backspace operation is done over the EOT, and 
then the current matrix heading, current matrix elements, and a 
new EOT is written. 
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A disk is preferred to a tape for the following reason. Be- 
cause of the physical separation of the read and write heads on 
most tape drives there may be tape tolerance problems thus back- 
spacing over the EOT is usually not successful. Instead of ending 
up positioned in front of the EOT, the write head is often posi- 
tioned in front of the previous matrix elements (E fl in the above 

sketch) . The current matrix heading will be written over the pre- 
vious matrix elements. This causes problems later when trying to 
read the records written on the tape. To alleviate this problem, 
it is strongly recommended that all FORMA tape subroutines (INTAP , 
LTAPE RTAPE, WTAPE , and UPDATE) use an intermediate device such 
as a disk. At the start of a computer run, the existing tape 
should be copied onto the disk by using computer control cards . 
Likewise, at the end of the run, the disk should be copied back 
onto tape by using computer control cards. 
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APPENDIX B - - BASELINE USER-PAK DEFINITION 


This appendix presents a listing of the seven subroutines and 
two functions which comprise the basic user-supplied packages 
required for any simulation. 
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